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1 .  INTRODUCTION 


We  present  the  results  of  an  extensive  investigation  of  the  electrical 
characteristics  of  the  concentric  spherical  cavity.  Although  this  system 
possesses  a  very  simple  geometry,  to  date,  an  analysis  of  this  cavity  ha3 
nevertheless  been  almost  totally  ignored.  The  only  meaningful  earlier  work  is 
a  short  note  by  Broc.1  Broc's  work  confined  itself  to  determining  the  Q  of  a 
silver-walled  concentric  cavity  at  a  fixed  frequency  for  each  of  two  low  order 
modes.  This  note  is  far  from  a  complete  solution  of  the  concentric  spherical 
cavity  problem.  Since  this  cavity  serves  as  the  starting  point  for  the  solu¬ 
tion  of  a  number  of  important  practical  electromagnetic  problems,  especially 
when  considered  in  one  asymptotic  limit  or  another,  it  was  necessary  to  care¬ 
fully  examine  it  in  its  simplest  form.  The  results  obtained  in  the  initial 
stages  of  this  study  are  rather  interesting  in  themselves,  since  they  give 
insight  into  electrical  characteristics  of  many  problems  of  a  considerably 
more  complex  nature.  Indeed,  a  series  of  subsequent  reports  will  elaborate  on 
a  number  of  such  important  practical  problems,  ranging  from  coupling  phenomena 
to  antenna  properties. 

At  the  outset  a  bookkeeping  problem  presents  itself.  To  deal  with  this 
situation,  the  results  of  our  investigations  will  be  discussed  in  a  number  of 
reports.  Each  report  will  be  sufficiently  restricted  both  in  size  and  subject 
to  render  the  whole  manageable  while  at  the  same  time  retaining  coherent 
meaning  and  connectivity  of  one  to  the  others. 

This  particular  report  in  the  sequence.  Part  I,  considers  the  geometric 
characteristics  of  the  concentric  spherical  system.  The  formal  vector  field 
solution  appropriate  for  the  system  appears  in  Part  II.  Part  III  contains  the 
formal  analysis  for  the  eigenvalue  solution  for  this  boundary  value  problem, 
and  Part  IV  gives  the  numerical  solutions  obtained  for  the  eigenvalues  and 
plots  of  their  trajectories,  along  with  a  discussion  of  the  results  presented. 


2.  FORMAL  SOLUTION  FOR  FIELDS  WITHIN  CAVITY 
2. 1  Basic  Field  Theory 

We  shall  assume  harmonic  time  dependence  of  the  form  exp(j<ot).  Hie 
explicit  fields  we  seek  are  those  appropriate  to  the  concentric  spherical 
cavity  illustrated  in  figure  1.  At  the  start  we  must  establish  some  conven¬ 
tional  notation;  hence,  we  denote  the  radius  of  the  outer  sphere  by  r  ■  a  and 
the  radius  of  the  inner  sphere  by  r  ■  b.  We  shall  find  it  convenient  to 
introduce  a  parameter,  namely,  the  ratio  of  the  inner  to  outer  sphere  radii, 
which  we  denote  by  R  *  b/a.  We  shall  restrict  our  considerations  in  this 
report  to  the  situation  where  there  are  no  sources  present  within  the  cavity 
(except  of  course  on  the  boundary  walls)  and  to  perfectly  conducting  walls. 
Also,  we  shall  consider  at  this  time  only  an  isotropic,  homogeneous,  linear 
dielectric  medium  filling  the  cavity.  At  first  we  assume  a  non-lossy 
dielectric.  In  section  3,  we  shall  extend  our  results  to  lossy  dielectrics. 

J.  Broc,  Electromagnet  is me,  Academie  des  Sciences,  230,  Coaptes  Rendus 
(January  1950),  199-9. 
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Figure  1 .  Geometry  of  concentric  spherical  cavity. 


In  our  present  problem  the  Maxwell  equations  are 


$  x  E  =  -jupH  ,  $  •  H  *  0 

$  x  H  =  J  +  ju>£ E  ,  $  •  E  =  p/e 

Since  H  is  solenoidal,  we  have  then  a  vector  potential  X  such  that 

*(f)  -  $  x  X(*)  , 


which  then  yields 


=  -;}u>yX(?)  -  $*(?)  f  ( 

where  is  a  scalar  function.  X  and  ♦( r)  are  connected  via  the  relation 

$  x  $  x  X  -  w2peX  ■  5  -  ju>e$*  ,  ( 


which  we  can  write  alternatively. 


V2X  +  u^peX  ■  -3  +  •  X)  +  ja>e$$ 


We  can  without  loss  of  generality  take  the  gauge  condition 


$  •  X  ■  -jue#  , 


which  then  gives  us  the  relation  between  the  vector  potential  A  and  the  cur¬ 
rent  sources  J,  if  there  are  any: 

(V2  +  0)2ye)A  =  -3  .  (7) 

He  may  consider  this  vector  potential  X(r)  as  the  generator  of  an  electric  and 
a  magnetic  field  given  by 


i,(f)  -  -:«■*(?)  *  ?[?  •  *(?)]  <8.) 

6,(f)  -  1  *  X(f)  (8b) 


Now  in  a  source-free  region, 
cavity,  we  have  $  •  E  *  0;  hence,  we 
fields. 


such  as  within  our  concentric  spherical 
can  have,  independent  of  the  previous 


E  =  -$  x  F  ,  (9) 

where  the  vector  potential  F  gives  rise  to  a  solenoidal  electric  field.  We 
then  obtain  from  the  Maxwell  equations  a  corresponding  magnetic  field, 

fi(?)  =  -jo.e*(?)  -  ,  (10) 

where  is  a  scalar  potential.  $  and  f  are  related  through 

$  x  $  x  |  -  u)2ne$  *  >  (11) 

or  alternatively  through 

(V2  +  w2jie)f  *  $($  •  I)  +  jwpfa  .  (12) 

Again  without  loss  of  generality  we  can  fix  the  gauge  by 

$  •  f  *  ,  (13) 

and  obtain  for  the  vector  potential  ? 

(V2  +  w2^)?  =  0  .  (14) 

ttiis  vector  potential  generates  an  electric  and  a  magnetic  field  within  our 
source -free  cavity  given  by 

l2(£)  ■-$*$(?)  (15a) 


M?)  -  1(1)] 


(15b) 


Since  the  fields  generated  by  A  and  by  F  are  independent,  we  have  as  the  total 
fields  in  our  source-free  cavity  the  sum  of  the  two  contributions 


l(r)  =  -$  x  *(£)  , 


and 


*fr)  -*“*(*)  ♦  =  *  **  **  ‘ 


(16a) 


(16b) 


Now  we  can  always  decompose  a  vector  field  into  a  component  trans¬ 
verse  to  the  radial  unit  vector  $r  and  a  component  along  $r.  He  shall  refer 
to  the  latter  as  the  longitudinal  component.  Thus,  the  electric  field  at  a 
point  within  our  cavity  can  be  formally  decomposed  in  this  fashion  as  follows: 

E(E)  .  Sr[Sr  •  1(f)]  -  Sr  x  [*,  x  1(f)]  .  (17) 

He  shall  now  proceed  to  construct  the  total  electric  and  the  total 
magnetic  fields  within  the  cavity  as  superpositions  of  just  such  transverse 
components  and  longitudinal  components.  This  requires  some  arithmetic. 
Consider  the  specific  vector  functions 

a( r)  =  ejA^r)  ,  !(r)  s  erFr(r)  .  (18) 

Note  that  it  should  not  be  expected  that  in  general  we  can  satisfy  the  rela¬ 
tions 


Sr  .  (V*A)  -  v2(Sr  .  A)  or  •  (V*F)  =  7*(Sr  •  F)  ; 


that  is.  Ay  or  Fr  are  solutions  of  the  scalar  Helmholtz  equation.  Hhat  we  do 
have,  in  source -free  regions,  is 

$x$x£-<upe£*  -jioe$$  ,  (19a) 

and 

$  x  $  x  ?  -  U)2uef  =  —  .  (19b) 

From  these  we  can  find  the  differential  equations  that  Ar(?)  and  Fr(?)  must 
satisfy.  Let  us  work  first  with  £  =  £rAr.  Now  some  arithmetic  exercises  are 
necessary. 


In  spherical  coordinates  the  gradient  operator,  we  recall,  is 


r  3r  r  39  r  sin  9  3r 


(20) 
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V.v.v. 


1  V,  /. 
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Then  from  the  following  identities  connecting  the  orthogonal  unit  vectors  in 
cartesian  coordinates  and  spherical  coordinates. 


we  find 


®x  =  ®r  s*n  ®  cos  ♦  +  ®8  cos  0  cos  $  -  e^  sin  $  , 

%y  *  %r  sin  0  sin  $  +  eg  cos  0  sin  $  +  cos  $  , 

®z  *  ®r  008  9  -  e©  sin  0  , 

®r  "  ®x  sin  9  cos  $  +  ey  sin  0  sin  $  +  ez  cos  0  , 

£q  *  £x  cos  0  cos  ^  +  Sy  cos  0  sin  $  -  %z  sin  0  , 

-*•  ♦  .  ♦ 

e^  »  -  ex  sin  $  +  ey  cos  $  , 
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(21a) 


(21b) 


e^  sin  0  , 

(22a) 

cos  0  , 

(22b) 

-(er  sin  0  +  ee  cos  0)  . 

(22c) 

In  turn  we  find 


r2  sin  0  99  "®  r2  sin2  0  3$ 


r  sin  0  3$ 


V 


We  thus  have  for  our  choice  of  £(r)  substituted  into  equation  (19a);  that  is, 
$  x  $  x  (erAr)  -  w2pe(erAr)  =  , 

or  the  more  explicit  relation,  with  k2  =  oj2ug 


-rfe  [nr?  Is sln  9  W  *  — h  J? +  k2K(i))  riire  I?) 

•  ^  fc  ♦  *  £  w  ♦  r=b  lM 

Then  if  we  choose,  as  we  are  free  to. 


, r +1  ^Ar(r) 

-3“e*(rJ  =  - 37- 


this,  indeed,  follows  from  equation  (6);  we  will  obtain  as  the  differential 


equation  for  Ar, 

■*  32ftr  1 

.  3_  / 

f  .  3Ar\ 

1  32Ar  +  k2Ar] 

r  3r2  r2  sin  0 

30  \ 

30  ) 

+  r2 

sin2  0  3<j>2  J 

With  the  observation 

that 

1_ 

r 

3^_ 

3r2 

[■£01  ■ 

~( 
3r-  \ 

“r )  *  r  1?  (“r)  ' 

we  can  conveniently  rewrite  equation  (28)  as 


r  0  Ar ( ) 

( V2  +  k2 )  =  0  . 


From  equations  (24)  and  (16b),  we  then  conclude  that  Ar,  the  solution  of^the 
differential  equation  (30),  will  generate  a  magnetic  field  transverse  to  er. 
This  will  be  referred  to  as  the  TM  field  solution  in  the  following  discus¬ 
sions.  In  precisely  the  same  manner,  we  can  process  ?(?)  of  equation  (18)  and 
obtain 


(V2  +  k2)  Fr( r)/r  =  0 


as  the  defining  differential  equation  of  the  scalar  function  Fr(£).  Also  we 
have  as  in  equation  (24), 


*  e0  3Fr  e7  3Pr 

7  x  (erFrJ  *  r  8in  e  "  r  ae 

which  is  transverse  to  £r. 


•  »  .  •  k  •  .  •  .V  .  •  ■  ’»  .  •  .%  .  » 


*  *  *  “  •>  *  •**•»•«*  *  **  b.  “  h  w  •  •  »  '  a  w  1,  »  1  *  k  *  .  *  W  *  .  *  V  •  a . 


Thus,  Fr  will  generate  an  electric  field  within  the  cavity  which  is 
transverse  to  £r  and  which  will  be  referred  to  as  the  TE  field  solution  in 
what  follows. 

Once  again  we  change  notation  by  defining  two  scalar  functions  $(r) 
and  f|*(  r) ,  such  that 

A(r)  =  r(Ar/r)  =  ri|>(r)  (33) 

and 

$(*)  =  r(Fr/r)  =  rij>(r)  .  (34) 


The  electric  and  magnetic  fields  in  terms  of  these  scalar  functions  are 


i(i)  -  -?  x  [?♦(?)]  ♦  x  ?  x  [;*(?)]  us) 

and 

*(?)  -i «  [H(J)J  .3^?  x  ? ,  [;♦(?)]  ,  ,36) 

where  the  TM  part  of  the  fields  is  generated  by  the  scalar  function  ^(r)  and 
the  TE  part  separately  by  $(r],  and  of  course  we  have 

(V2  -  Y2)  ♦(*)  -  0  *  (37) 

and 

(V2  -  Y2)  ♦(*)  =  0  ,  (38) 

and  where  we  choose  to  replace  k  by  jy. 


It  is  useful  to  write  out  explicitly  the  field  components  in  terms 
of  the  scalar  functions  and  $(r).  We  have  then  from  equations  (35)  and 
(36): 


*(*)  - 

*»({)  ■  ~7TiT?ly  +  ’ 

*(*>  *  7  h  *  377  r~SInT  7H7  M*)]  • 

*(*)  -  3=7  [fe  -  ’2]  . 

"«(f)  '  3ST  r  9756  ♦  7-^TT  I?  • 

*(*)  *  377  7"SIiT7  M*)  J  -  7  If  ('♦<*'))  • 
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►V 


m 

1  / 
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2.2  The  Eigenmode 3 


As  usual  we  take  the  scalar  functions  as  products  of  radial  and 
angular  function  factors;  that  is, 

i}i(r)  -  i|>(r)Y(6,$)  (45a) 


and 


♦(*)  *  <fr(r)Y(9,<M  .  (45b) 

From  the  forms  of  equations  (37),  (38),  and  (45),  we  find  that  we 

shall  have  solutions  of  the  form 

'J'nmH*)  “  <'n<r>pnm(cos  0>  008  “♦  '  (46) 

’J'nm^r)  =  ’Mr>pnm(co8  0)  8in  * 


and 

♦nSH*)  *  <Mr>pnm(cos  9 )  cos  m$  ,  (47) 

♦nSH^)  “  *n<r>pnm<c°8  0>  sin  * 

where  n  is  a  positive  integer,  m  can  run  from  0  to  n,  and  the  functions 
are  the  familiar  associated  Legendre  functions  of  the  first  kind. 

The  radial  factors  satisfy  the  differential  equations 


[r*n<*>]  -  [f2  ♦  n(nr2  n]  l**n<*>]  “  0  (48) 

and 

[r*„(t)l  -  [r2  *  [r,„(D]  -  0  .  (49) 

The  formal  fundamental  solutions  of  equations  (48)  and  (49)  are  the  modified 
spherical  Bessel  functions  in(yr)  and  kn(yr),  where  the  in  are  singular  at 
r  «  *  and  the  kn  are  singular  at  the  origin.  Complete  information  concerning 
the  properties  of  these  functions  is  available  in  the  mathematics  reference 
literature.2 


In  the  concentric  spherical  cavity,  the  region  of  interest  excludes 
the  origin;  hence,  the  radial  factor  must  be  a  linear  combination  of  the  in 
and  kn  functions,  e.g. , 


2A.  Erdelyl,  ed .,  Ch  VII,  Higher  Transcendental  Functions,  California  Insti¬ 
tute  of  Technology,  Bateman  Manuscript  Project,  II,  McGraw-Hill  Book  Co., 
Inc.,  NY  (1953). 
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In  the  above  we  have  ma£e  another  attempt  to  condense  the  notation  by  defining 
for  our  vector  field,  V'J*  (r) ,  the  convention  that  j  ■  1  implies  TE  modes 
and  j  =  2  1M  modes. 

The  vector  field  eigenmodes  given  in  equations  (54)  and  (55)  form  a 
complete  orthogonal  set  for  our  cavity.  The  orthogonality  is  shown  in  appen¬ 
dix  A.  Any  electric  field  within  the  cavity  can  be  expanded  in  this  set  of 
eigenvectors.  Thus,  we  may  write 

t(J)  -  I,  I,  L  1  (»&°’S&0>(?)  ♦  .  (58, 
n=1  p*1  m=0  j  =1  '  > 

where,  integrating  over  the  cavity  volume. 


where 


4^,S)  =  ///  *(*)  •  ^i8)(t)  *  [4U8)]2  ■ 

k 

[flL|i§)]2  8  ///  . 


(59a) 


(59b) 


We  have,  of  course,  the  vector  relations 


? « ? « iurh f)  - 

*  * « « 4&$>m  * 


and  we  have,  in  addition,  the  relations  between  the  electric  and  magnetic 
eigenvectors 


?  *■  ■*£&%)  • 


Any  magnetic  field  within  the  cavity  similar  to  equation  (58)  can  be 
written  as  an  expansion  in  the  complete  set  of  magnetic  field  eigenmodes  of 
equations  (56)  and  (57).  This  in  its  turn  can  be  written  using  equation  (62) 


Lilli-2* 

n-1  p-1  m-0  j-1  \ 


pm  j 

0TT2 


=  I6(r  -  r') 


I  I  I  I 

n-1  p-1  m-0  j »' 


where  I  is  the  unit  dyad.  We  can  show  that 


^njlm  ^  =  ^n^m  ^  • 


and  we  therefore  need  consider  only  one  set  of  normalization  constants,  which 
we  shall  assume  to  be  the  ^npm°  *  The  proof  is  straightforward  using  the 
vector  identity 


?  •  (M  X  5)  =  li  .  ?  x  5  .  s  .  ^  x  s 


and  choosing 


*  .  >(*)  ,  hhlij;  >(*)  . 


(68a) 


(68b) 


We  get  upon  substitution  into  the  identity  and  integrating  over  the  volume  of 
the  concentric  spherical  cavity 


HI  dT(f)  ?  •  (*&  >(*)  «  «  «  *&  >(*)] 

-  ///  ®(t) » .  ijjs,  >(t)  •  $  *  t<k  >(*) 

-  ///  4t(t)  »(*)  •  ♦  * » •  »(*)]  • 

Rearranging  and  applying  Stokes  theorem  this  becomes 

in  dt(?)  ? » ’(i) .  t « * ^  \t) .  //vjj(?)  •  fep:  ’({) 

walls 

>  * » >(*}}  -  frij’]2 111  d,(*)  •(*)  •  »(*)  • 


(68c) 


( 68d) 


The  surface  integral  vanished  due  to  the  vanishing  of  Etan  on  the  cavity 
walls.  Then,  by  equation  (62),  we  can  rewrite  this  as 


-h&’]2  ///  «(?)  ’C?)  •  ’(i) 

-  -[^’]2  ///  «(t)  I &  '(?)  •  '(f)  - 


(68e) 


from  which  the  identity  in  equation  (67)  follows. 

Weyl3  demonstrated  earlier  that  a  dyadic  Green's  function  exists  for 
the  geometry  we  are  dealing  with.  This  function  is  of  considerable  signifi¬ 
cance.  For  the  electric  field,  the  dyadic  Green's  function  Gg(r,r')f  which 
satisfies 


?  X  ^  X  GE(r,r')  +  Y2GE(r,r') 


(69a) 


is  given  by 


»npm 


Jnp 


(69b) 


Weyl,  Uber  die  Randwertaufgabe  der  Strahl^ungstheorie  and  Asyaptotische 
Spektralgesetze,  Z.  Reine  u nd  Angewand .  Nath,  143  (1913),  177-202 . 
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'  w  W  -  A 


rtrr*: 


' .V TTTS 7  ’>  ■> ■I*?  '5 "A  ’.*  W "X ■>  ■.« ■« 

,'>  ,'"V 
•  -.  »\. 


As  a  demonstration  we  show  that 


ft 


///  dx(r')  E(r')  •  ($♦,  x  $+,*  GE(r,r')  +  y2GE(r,r')  +  I«(r  -  r')j 

r*r\r»r*  _  cnhpr  _  *  ^  * 


cone,  spher 
cavity 


-iii  I  (^'J^C///  ar(io  S(?.) .  sA|ie,(?o]  ^e>c?)[oAJie)] 

n=1  p=1  m=0  j=1  * L  J 

*  [///  d,(f.)  S(f')  •  ^i0>(?)[oifc0,]-2|  *  [r2  -  (rij>)2] 

-t2  I  i  i  l  ([///  «(M  «(?•)  • 

n*1  p=1  m=0  j=1  ‘ <- 

♦[///  dt(r')  S(f')  •  ^°>(f,)]'iji0,(J')[aiU0>]'2)  *  -  (''ip’)2] 

J  \  D- 


-2 


.///  dx(f')  £(f')d(f  -*•) 
-  -E(f)  +  E(J)  -  0  . 


We  note  that  the  poles  of  the  Green's  function  for  the  electric 
field  given  in  equation  (69)  are  the  eigenvalues  or  resonances  of  our  perfect 
cavity. 

The  magnetic  field  dyadic  Green's  function  is  given  by 

1-2 


-  -  l  l  If 

n=1  p»1  m»0  j»1  '  J 


♦  ^■”(?)Si|;0,(*r-)[A^0>]'2)  ♦  [r2  -  (ril1)2] 

which  satisfies 

$  *  $  *  G|j(r,r ')  +  Y2SH(r,r')  -  -lfi(r  -  r')  . 


(70) 


(71) 


Having  formally  discussed  the  vector  fields  for  the  cavity,  we  proceed  next  to 
a  formal  discussion  of  the  eigenvalues. 
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2.3  The  Eigenvalues 


He  shall  discus  only  briefly  the  determination  of  the  eigenvalues 
Ynji^  which  yield  the  eigenfunctions  ^np(r)  and  4>np(r)»  It  is  somewhat  sim¬ 
pler  to  consider  the  TE  fields  which  involve  the  scalar  functions  $np(r).  At 
the  cavity  walls,  the  boundary  conditions  that  the  electric  fields  are  requir¬ 
ed  to  satisfy  are 

E0;nm<r  =  a)  =  0  =  E0;nm(r  =  b)  ,  (72a) 

and 

E$;nm(r  =  a ^  ~  ®  =  E^;nm^r  ~  b)  •  (72b) 

Equation  (40)  then  requires  that  we  have 

4>n(a)  =  0  ,  (73a). 

and 

Vb)  =  0  '  (73b) 

which  can  be  satisfied  if  we  take  4>n(r)  in  the  form 

♦n(r)  -  in(Yr)kn(Yb)  -  in(Yb)kn(Yr)  .  (74) 

Clearly,  equation  (73b)  is  automatically  satisfied.  To  satisfy  equation 
(73a),  we  must  find  the  roots  or,  equivalently,  those  values  of  y  for  which 

in(Ya)kn(Yb)  -  in(Yb)kn(Ya)  -  0  .  (75) 

In  section  2.4  we  shall  go  into  great  detail  concerning  the  analytic 
procedure  for  solving  this  eigenvalue  problem.  He  shall  see  at  that  time  that 
there  is  a  discrete,  ordered  infinite  set  of  roots  for  each  order  n.  These 
roots,  which  we  denote  by  y  *  jYnp^  f°r  P  =  1  •  2,  ...»  lie  along  the  positive 
imaginary  axis  in  the  complex  y -plane.  For  each  of  these  values  we  obtain  a 
radial  factor;  i.e., 


♦np<r>  2  infjYip’rJjCn^Yip’b)  -  infjyip^j^n^Tip’r) 


The  functions  given  in  equation  (76)  used  in  equation  (52)  will 
generate  the  corresponding  TE  vector  eigenmode  fields  from  equations  (40) 
through  ( 44 ) . 

Por  the  TM  modes,  we  need  to  determine  the  eigenfunctions  <|;  (r) 
which  satisfy  the  boundary  conditions  given  in  equation  (72).  From  equations 
(41)  and  (46)  we  find  we  must  have 


MXl-V.  VATV  w  vTOTvVT*. 


It  i 


^nD<r). 


at  r  *  a 


Is-  lrtnp(r)]J  -  0 


(77a) 


(77b) 


at  r  *  b  . 


These  requirements  can  be  met  if  we  take  ^n(r)  of  the  form 

♦n(r)  -  in(Yr)[Ybkn(Yb)]'  -  kn(Yr)[Ybin(Yb)]'  ,  (78) 

where  [XZn(X)]'  means  differentiate  with  respect  to  the  variable  X.  We 
observe  that  equation  (77a)  is  automatically  satisfied.  Equation  (77a)  re¬ 
quires  that  we  have 

[Yain(Ya)]'[Ybkn(Yb)]'  -  [Yakn(Ya)]'[Ybin(Yb)j*  -  0  .  (79) 

We  shall  see  in  the  following  section  that  this  relation  will  be  satisfied 
only  for  an  ordered  infinite  discrete  set  of  values  y  -  Y^n  ’  for  p  *  1,  2, 
. . . ,  for  each  order  n.  The  physically  meaningful  roots  lie  along  the  positive 
imaginary  axis,  just  as  for  the  TE  case.  Each  such  eigenvalue  Y^p  gives  rise 
to  a  radial  factor  function,  i.e.. 


♦npIO  5  in(jrjp’r)  [ir^’b  k„(jY ‘£’b)] 

-  ^^[r^’b  ln(<'b)]'  • 


Functions  given  in  equation  (80)  used  in  equations  (50)  generate  the  corre¬ 
sponding  TM  vector  eigenmode  fields  from  equations  (39)  through  (41)  and  (43) 
and  (44). 

Next,  we  briefly  discuss  the  effect  of  a  lossy  dielectric  on  these 
eigenvalues  and  on  the  corresponding  vector  eigenmodes. 

2.4  Effect  of  Lossy  Dielectric  on  Eigenvalues  and  Eigenmodes 


By  allowing  the  dielectric  to  have  a  finite,  isotropic,  homogeneous, 
frequency  independent,  linear  conductivity,  we  can  account  for  the  effects  of 
a  lossy  dielectric  filling  the  concentric  spherical  cavity.  We  now  have 


5(?,t)  *  ol(r,t)  . 
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k a 


,  V  J"  ST-.  •-.  r. 


The  Faraday's  Law  and  Ampere's  Law  equations  can  now  be  written 


♦  +,+  .  3H(r,t) 

?  x  E(r,t)  .  -a  — JJ-  . 


$  x  H(r,t)  =  oE(r,t)  +  e 


from  which  we  obtain  the  relation 


3*(r,t) 


3t 


(82a) 

(82b) 


&Cr.  t)  attf.t)  ***,*, 

ye  - - —  +  0y  — - -  +  vr  X  7r  X  E[r,tJ  a  0  . 


at* 


(83) 


Since  we  have  the  same  boundary  conditions  for  £(?,t)  to  satisfy ,  we 
can  expand  it  in  the  same  complete  orthogonal  set  of  electric  eigenvectors 
given  in  equations  (54)  and  (55).  Equation  (83),  then,  is  equivalent  to 


I  l  l  I  [ue  ///  dx (?•) 

n=1  p=1  m=0  3=1 


at2 


jses. 


[0&e’]2  [»&3,12 


*  TO  III  at (?•) 


.  ,  U!  atfr  v 

I  [4k'\2  [s&’l2  I 

.  *&0,(t-)*r»M&.0,(i)  „ 

[aili*’]2  [oi|i0,J2 

If  we  dot  through  with,  say,  and  integrate  with  respect  to  ? 

the  cavity,  this  reduces  to 

»*  £-  US  dx(r)  l(;,t)  •  4k  >(1) 

3t2 

-  TO  it  ///  at(i)  S(i,t) .  4k  ’(?) 

-  [T^ipl2  SSI  «(f)  *(?,t)  •  lili  >[?)  -  0  , 

where  yQ  is  the  zero-conductivity  value. 
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over 


(84) 


Assuming  time  dependence  of  the  form  eJw  ,  where 

to  *  <o1  +  ju)2  * 

we  can  solve  for  the  real  and  imaginary  parts  of  the  frequency, 
then  becomes 


Equation  (84) 


new2  -  jouu)  +  [Yo?np]2  =  0  ' 


where 


Yo^np  -  j^o.np*^2  » 


and  u0fnp  is  the  eigenf requency  for  a  «*  0: 


“np  =  “o,np 


“  (^o,npe)' 


Adopting  the  conventional  notation  for  the  quality  factor  Q  of  the 
system,  we  define  for  the  complex  eigenf requency,  when  o  *  0: 


„  _  “o,npe 

Qdiel, np  -  Z 


Then  we  rewrite  the  eigenfrequency  for  a  simple  lossy  dielectric 
filling  the  cavity  as 


“np  “  “o,np  71  “ 


( 22diel, np)2  +  3  22diel,np 


The  imaginary  part  of  the  eigenf requencies  introduces  a  constant  uniform 
damping  factor  over  the  entire  cavity  volume,  namely  exp(-ot/2e),  independent 
of  frequency.  This  same  damping  is  applied  to  every  vector  eigenmode.  As  the 
conductivity  of  the  dielectric  increases  from  zero  we  note  that  not  only  does 
the  damping  factor  become  stronger,  but  the  frequency  of  oscillation  decreases 
for  each  eigenmode  until  it  eventually  becomes  zero  at  a  *  2ew0  np  for  the  n, 
path  mode.  ' 

The  trajectory  can  be  quite  readily  established  for  the  eigenvalues 
<0  in  the  complex  plane  as  the  conductivity  increases  from  zero.  Multiply 
equation  (90)  by  j/“0fnp  to  get 


\“o,np/ 


j  /I  - 


(2Qdiel,np)2  "  22diel,np 
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Define 


x_  -  Mo>nP  =  £_ 
"  22diel, np  2e 

and  we  can  write  equation  (91)  as 


which  we  can  rearrange  to 


From  this  relation  we  conclude  that  each  eigenfrequency  starts  out  on 
the  positive  imaginary  axis  at  zero  conductivity.  As  the  conductivity  in¬ 
creases,  each  eigenfrequency  moves  along  a  quadrant  of  its  own  circle  in  the 
complex  plane  with  its  center  at  the  origin.  At  a  =  2Ea,o,np'  the  e*-9envalue 
reaches  the  negative  real  axis  and  all  oscillatory  temporal  behavior  ceases. 
All  fields  at  this  value  of  o  merely  damp  out  exponentially  with  time. 

Thus  far  we  have  considered  only  the  effect  of  a  simple  lossy  di¬ 
electric  on  the  cavity  modes.  In  a  subsequent  report  we  shall  extend  the 
discussion  to  include  the  more  realistic  case  of  finite  conductivity  for  the 
boundary  walls. 

We  now  proceed  to  a  detailed,  look  into  the  method  used  to  find  the 
eigenvalues  explicitly. 


3.  ANALYTIC  AND  NUMERICAL  DETERMINATION  OF  THE  EIGENVALUES 


General  Range  of  the  Eigenvalues 


We  shall  now  present  detailed  analyses  of  the  methods  used  to  deter¬ 
mine  the  numerical  values  of  the  eigenvalues  for  a  nonlossy  dielectric 
within  the  cavity.  The  effect  of  a  lossy  dielectric  has  already  been  dis¬ 
cussed.  The  analyses  to  be  presented  are  quite  tedious.  Nevertheless,  it  is 
significant  to  discuss  this  in  some  depth,  if  only  to  put  proper  perspective 
on  why  the  concentric  spherical  cavity  problem  has  only  now  been  resolved  in 


detail.  In  our  investigation  of  the  eigenvalue  problem,  we  evaluated  all 
those  y’s  such  that  Y^p  S  2 ®  f°r  n  *  *  through  4  in  the  TE  case  (j  *  1)  and 
for  n  *  1  through  5  in  the  TM  case  (j  ■  2).  For  the  nonlossy  dielectric  the 


physically  acceptable  solutions  for  the  eigenvalues  lie  along  the  positive 
imaginary  axis  in  the  complex  y  plane.  This  will  become  evident  as  we  proceed 
through  the  following  discussions.  The  eigenvalues  are  determined  as  func¬ 
tions  of  R,  i.e.,  the  ratio  of  inner  to  outer  radius. 


3.2  The  TE  Eigenvalues  yAp' 

It  is  simpler  to  begin  with  the  analysis  for  the  determination  of  the 
TE  eigenvalues.  We  indicated  earlier  that  for  j  =  1  we  need  to  find  the  roots 
for  the  functional  relation  in  equation  (75),  which  we  repeat  here: 

in(Ya)kn(yb)  -  in(yb)kn(ya)  =  0  , 

For  convenience,  we  introduce  the  notation 

jx  =  ya  ,  (93a) 

jy  =  Yb  =  yRa  =  Rjx  ,  (93b) 


and  hence  what  we  seek  are  the  values  of  x  for  each  value  of  n  such  that 

in( jx)kn( jRx)  -  in(jRx)kn( jx)  =  0  .  (94) 

Note  that  from  the  properties  of  the  modified  spherical  Bessel  functions, 
equation  (75)  will  only  possess  roots  along  the  imaginary  axis  in  complex 
conjugate  pairs.  Only  those  on  the  positive  half  of  this  axis  correspond  to 
physically  acceptable  solutions  for  our  cavity  modes. 

3.2.1  The  YU*  (R) 

For  n  =  1  the  explicit  modified  spherical  Bessel  functions  are,  for 

real  x. 


i,<ix)  ,  ai"h<3«>  .  .  J: 

(jx)*  Jx  ' 


sin  x  cos  x 


)  • 


,  .  .  »  -jx/  1  1  \  /cos  x  .  sin  a  .  ./sin  x  cos  x\ 

k’(3X)  *  •  fen?  *  <w)  ■  -(—  *  -sr^ +  — — ) 

Substituting  these  forms  into  equation  (94)  for  n  =  1  gives 


sin  x\ 

/sin  y 

cos  y' 

\  /sin  x  cos  x\ 

/cos  y  .  sin  y\i 

x  / 

V  y2 

y  , 

’  V  x2  *  / 

\  y2  •  Y  // 

(95) 

(96) 


0  . 


Note  that  the  real  part  of  the  expression  is  identically  zero.  After  some 
algebra,  we  can  rearrange  this  to 


(« *  57) ,in(x  - y)  ■  Cjr  •  -k) oo,(x  - yl  ■ 0  • 

or  even  more  conveniently  to 

tan [ ( 1  -  R)x]  -  (1  -  R)x/[1  +  Rx2 )  H  f1(R>x)  .  (97) 
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It  should  be  clear  that  our  problem  of  determining  the  eigenvalues 
here  reduces  to  finding  the  intersections,  for  given  R,  of  the  rational  frac¬ 
tion  f i (R;x)  with  the  tangent  function.  We  shall  see  that  this  will  indeed  be 
the  case  for  all  orders .  ,i\»  The  intersections  can  be  uniquely  ordered  in 


:<A‘, 


increasing  value,  i.e.,  xj'{(R)r  xj^2(R),  X1,3^R)# 


We  further  observe  that  since 


fi<R;x)  =  (1 

-  R) ( 1  -  Rx2)/(1  +  Rx2)2  , 

(98a) 

fj(R;0)  =  (1 

-  R)  , 

(98b) 

we  have  both  the  tangent  function  and  its  derivative,  respectively,  coincident 
with  fi  and  its  derivative  at  the  origin.  We  further  note  that  f1  is  never 
negative  and,  furthermore,  goes  to  zero  asymptotically  with  x  ♦  ®.  Hie  only 
zero  of  fj  occurs  at 


Xq  *  •  (99a) 

Since  f^  starts  at  the  origin,  goes  positive,  and  ultimately  goes 
to  zero  as  x  >  ®,  equation  (99a)  locates  the  one  and  only  maximum  for  f^, 
where  we  have  actually 


f i (r;xq)  -  1/2  /l-  -  /l\  -  1/2  -1-£-r)- 
\/r  /  /r 

Also,  we  readily  find  for  R2  >  R1 , 


(99b) 


1//57  >  1//R2 


(100a) 


fl(Ri;x)  ( 1  -  R-|)  ( 1  +  R2x2] 
fi(R2;xJ  "  (l  -  R2J  (i  +  riX2)  >  1  ' 

because  1  -  R2  <  1  -  R1  and  1  +  R2x2  >  1  +  R^x2.  Thus,  we  see  that  as  a 
function  of  R  for  fixed  x  we  have 

fl(R,;x)  >  f1(R2;x)  >  f1(R3;x)  >  ...  for  R,  <  R2  <  R3  ,  (100b) 

and  the  maximum  value  occurs  earlier  for  the  larger  value  of  R.  We  illustrate 
this  in  figure  2.  Note  (1  -  R)  is  the  initial  slope  at  the  origin.  We  now 
know  all  the  essentials  of  the  characteristics  of  the  family  of  functions 
fj(R;x).  Let  us  examine  the  family  tan(1  -  R)x  in  some  detail.  For  a  given  R 
this  function  has  singularities  at 
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Figure  3.  Family  of  tan(1  -  R)x  for  R,  <  R2. 


Hence,  for  x  just  greater  than  zero,  f.|(R;x)  lies  below  the  tan(1  - 
R)x  function.  This  is  the  relative  position  of  the  two  functions  as  x  in¬ 
creases  out  to  the  first  singularity  of  the  tangent  function.  Thus,  we  con¬ 
clude  that  there  will  be  no  intersection  of  tan(1  -  R)x  and  f 1 (R;x)  for  x  < 
ir/2(1  -  R);  l.e.,  before  the  tangent  attains  its  first  singularity.  Further¬ 
more,  the  tangent  function  returns  from  at 

n 

X  “  2(1  -  R) 

and  remains  negative  as  x  increases  until  x  *  ir/(1  -  R).  The  first  crossing 
of  f.j(R;x)  with  the  tangent  then  occurs  at  x  >  ir/(1  -  R).  Actually,  x{^  (R) 
will  occur  between  x  =  tt/(1  -  R)  and  x  =  3ir/2 ( 1  -  R),  where  the  tangent  is 
increasing  with  x  from  zero  to  +■>.  All  further  crossings  will  occur  on  the 
arcs  of  the  tangent  function  where  it  is  positive;  i.e.,  xjp^(R)  will  lie  in 
the  range  pn/(1  -  R)  <  x  <  (2p  +  1)ir/2(1  -  R)  for  p  =  1,  2,  ....  For  a  given 
ratio  R  as  p  increases,  the  crossings  are  approximately  periodic  with  the 
period  of  the  tangent  function.  The  larger  the  values  of  p,  the  closer  this 
periodicity  becomes.  Since  the  tangent  function  essentially  spreads  out  with 
higher  values  of  R,  the  successive  crossings  occur  at  larger  intervals.  We 
have  already  noted  that  f^(R;x)  flattens  out  as  R  increases.  Then  we  conclude 
that  the  crossings  will  occur  closer  and  closer  to  the  zeros  of  the  tangent 


function  with  increasing  R,  which  we  have  already  noted  are  periodic  in  x  with 
increasing  period  as  R  increases.  The  details  of  the  behavior  of  the  xjp'(R) 
are  rather  nicely  depicted  in  figures  4  through  12,  which  show  the  two  func¬ 
tions  for  x  out  to  20  and  for  0.1  ^  R  i  0.9  in  steps  of  0.1.  Note  that,  for  R 
»  0.9,  xj]'(0.9)  occurs  beyond 


1  -  R 


10ir  >  20 


and  also  the  tangent  attains  its  first  singular  value  at 


2(1  -  R) 


Sir  =  15.7080 


(see  fig.  12),  which  is  the  only  singularity  or  zero  in  the  range  of  x  to 
which  we  have  restricted  our  interest. 
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12.  Generating  functions  for 
with  R  ■  0.90. 
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3.2.2  The  yU'(R) 


The  modified  spherical  Bessel  functions  for  n  =  2  are 


=  r 3  +  m 

sinhHir)  3  cosh<3x) 

T,  cos  x  3 

1_(  j X  ) 3  3XJ 

(jx)2 

L  x2  x3 

=  e-jxr  3  + 

3  +  ’1  -  [V3  -  M 

...  ..  3  cos  x| 

"  L(jx)3 

(jx)2  3XJ  [_\x3  x/ 

sxn  x  _  I 

J 

(102) 


.  |~3  sin  x  .  /3  1  \  1 

3  hi?-  *  fer  - ;)  008  *|  • 


(103) 


which  give  for  n  =  2  in  place  of  equation  (94),  upon  substitution  and  simpli¬ 
fication. 


° '  •(?'«)  Sl"  *][?  8l"  y  +  '  y)  003  y] 

'  -(?-?) sin  y][? 8l"  *  *  (?  -  i) 008  *]}  • 


(103a) 


We  note  again  that  the  real  part  of  the  expression  is  identically  zero.  With 
some  algebraic  exercising,  this  relation  can  be  put  in  the  form 


[(?  -  -  y)  *  8ln(x  ‘  y) 

'  “  y)  '  ^ "  *)] cos<x " y> 


or,  in  turn,  into  the  more  convenient  form 


tan[(1  -  R)x] 


_ (1  -  R)x(l  +  -j  x 

1  -  (R2  -  3R  +  Of2  + 


(103b) 


=  f2<R;x) 


(104) 


Again,  as  for  the  previous  case  for  n  *  1 ,  we  have  reduced  our 
eigenvalue  problem  to  determining  the  values  of  x^p^(R),  which  now  correspond 
to  the  intersections  of  the  tangent  function  on  the  left-hand  side  of  equation 
(104)  with  the  rational  fraction  function  f2(R;x).  Consider  the  denominator 
of  f2(R;x).  We  immediately  see  that 
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V 

k . 
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S'. 
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(a)  For  R  >  RQ  =  ( 3  -  / 5)/2  =  0.381/966,  the  denominator  is  posi¬ 
tive  for  x  >  0.  In  turn,  f2(R;x)  will  be  non-negative;  starting  at  zero  at 
the  origin,  it  will  at  first  increase  with  x.  As  x  ♦  •,  fj  +  0  asymptoti¬ 
cally. 


(b)  For  R  =*  Rq  the  denominator  has  a  double  real  root  at 

x(Rq)  =  x0  =  ,/3/R0  =  2.802,517  .  (105) 

It  then  follows  that  f2(Ro;x)  is  non"n®9®tive  starting  at  zero  at  x  *  0, 
becomes  infinite  at  xQ,  and  also  goes  to  zero  asymptotically  as  x 

(c)  For  R  <  Rq,  the  denominator  has  two  distinct  real  positive 

roots  at 

x_(R)  =  /3(r2  -  3R  +  l)/2R2  (l  -  /l  -  [ 2R/(r2  -  3R  +  l)]2}1^2  *  (106a) 

x+(R)  »  / 3(R2  -  3R  +  1 )/2R2  |l  +  /l  -  [ 2R/(R2  -  3R  +  l)]2Jl/2  ,  (106b) 

where  we  should  note  that  x_(R)  <  xQ  <  x+(R)  and  that  as  R  ♦  RQ 

x+(R)  -  x_(R)  =  2  / 3( R2  -  3R  +  l)/2R2  A  -  [2R/(R2  -  3R  +  l)J2  +  0  (107) 

because  the  right-most  factor  vanishes  in  this  limit.  The  function  f2(R;x) 
starts  out  as  zero  at  x  =  0,  becomes  positive,  and  increases  as  x  increases  at 
first.  f2  •  as  x  ♦  x_.  It  is  negative  between  x_  and  x+.  As  x  increases 
from  x_,  f2  returns  from  then  again  goes  to  -»  as  x  ♦  x+.  For  x  >  x+,  f2 
returns  from  +»  to  eventually  go  asymptotically  to  zero  as  x  -►  ».  For  x  >  x+, 
f2>  0. 


One  further  consideration  is  the  case  for  R  «  1.  Here  we  find  x_ 
"  /T  and  x+  *  vTT/R  »  1. 

We  have  thus  located  all  the  zeros  and  singularities  of  f2(R;x). 
Additional  information  can  be  obtained  from  f2(R;x),  which  is 


—  ~R-  *2(R;x> 


{l  +  (R2  +  l)(x2/3)  -  R(R2  +  l)(x2/3)2  -  R3(x2/3)3} 
{l  -  (R2  -  3R  +  l)(x2/3)  +  R2(x2/3)2}2 


(108) 


Note  that  the  denominator  is  never  negative.  For  case  (a),  the  denominator 
has  no  real  roots.  Upon  simple  test,  the  numerator  which  is  a  cubic  in  x2 
reveals  that  it  can  have  at  most  one  real  positive  root  for  x2.  Since  the 
function  f2  initially  has  a  positive  slope  near  x  ■  0  and  a  negative  slope  for 
large  x,  we  conclude  that  there  is  one  and  only  one  real  positive  root  for 
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f$(R?x)  in  case  (a).  Thus,  for  case  (a),  i.e.,  R  >  rq,  we  see  that  f2(R?x) 
behaves  essentially  like  f^(R;x),  which  we  examined  earlier.  Now  in  case  (b). 


the  denominator  of  f,  has  a  double  root  at  x 


numerator  in  equation  (108).  It  also  has  only  a  single 
=  / 3/R  ,  and  hence  we  can  write 


Consider  then  the 
root  which  occurs  at  x 


1  -  R. 


f2(Ro?x)  = 


{'  *  K  ♦  ♦  0^2  -  »|^j 


fx2  1  I* 

r'^J 


(*o)‘ 


We  see  that  f^(RQ;x]  is  singular  at  x  =  yS/Rq.  Also  for  x  <  xQ,  f^(RQ;x)  >  0 
and  for  x  >  xQ,  f^[RQ;x]  <  0.  Since  f2(Ro'"x)  has  no  zeros  for  x  *  xQ,  it  is 
monotonic  on  either  side  of  xQ.  Next  we  must  give  further  effort  to  examine 
f%2  for  case  (c),  where  R  <  RQ.  Again  the  numerator  in  equation  (108)  has  at 
most  a  single  real  positive  root.  This  must  occur  somewhere  in  the  range 
x  <  x  <  x+,  where  f2  has  to  change  from  positive  to  negative  to  permit  f2  to 
behave  properly  there.  In  figure  13  we  illustrate  f2(R;x)  for  R  <  RQ. 


Figure  13.  Rational  fraction  function  f2(RQ;x) 
for  R.j  K  R2  K  Rq* 
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Next  we  consider  the  crossings  of  f2(R;x)  and  tan(1  -  R)x.  In  the 
neighborhood  of  the  origin ,  the  terms  to  order  x®  of  the  series  expansions  of 
these  functions  are  given  by 

f 2(R;x)  *  (1  -  R)x  +  (1  -  R)3  ~  +  5{(1  -  R)b  -  R(1  -  R)[R2  +  (1  -  R)]}  ^  , 


tan( 1  -  R)x  =  (1  -  R)x  +  (1  -  R)3  j-  +  ^  (1  -  R)5x5  . 

Subtracting  coefficients  of  the  x^  terms,  we  find  for  R  <  1  that  the  differ¬ 
ence  is 


(1  -  R)5  -  {(1  -  R)5  -  R(  1  -  R)[R2  +  (1  -  R)]} 

-  — +  i  R<’  -  r)[r2  +  n  -  R)]  >  o  . 

We  conclude  then  that  f2  lies  below  the  tangent  function  in  the  neighborhood 
of  x  =*  0.  For  case  (a),  R  >  RQ,  the  story  is  essentially  that  for  the  TE,  n  * 
1  case.  In  case  (b),  where  R  =  R_  we  have 

(2(Ro,„)  .(,  -  «j.(.  -  r  "2)2  • 


The  first  singularity  of  the  tangent  is  at  x 


2.541,602,  which  occurs  before 


2.802,517.  Hie  tangent 


w/2(l  -  Rj 

igent,  after 


increases  beyond  this  first  singularity,  is  negative  out  to  it/(l  -  RQ)  ■ 
5.083,203.  Thereafter,  it  becomes  positive  and  increases  with  x.  Then  in 
case  (b),  the  first  crossing  occurs  after  x  ®  */(l  -  RQ),  where  f2(R;x)  is 
decreasing  from  and  the  tangent  is  increasing  from  zero  to  its  next  singu¬ 
lar  value  at  3ir/2(l  -  RQ) .  This  is  illustrated  in  figure  14.  The  successive 
crossings  will  occur  on  the  successive  positive -going  arcs  of  the  tangent 
curve.  The  higher  order  ones,  x^^Rq),  will  occur  at  close  to  the  period  of 
the  tangent  curve. 


Finally,  we  have  a  look  at  case  (c),  where  R  <  RQ.  Again,  f2(R;x) 
and  tan(1  -  R)x  coincide  and  have  the  same  derivatives  at  the  origin.  In  the 
neighborhood  of  x  *  0,  we  showed  that  tan  (1  -  R)x  sits  above  f2(R;x)  and  both 
functions  are  positive  and  increasing  monotonically  until  they  become  singular 
at  their  respective  singular  points.  Let  us  compare  the  series  expansions  of 
the  two  functions: 


-  -~R-  — ■  tan(1  -  R)x  *  1  +  (1  -  R)2x2  +  -|  (1  -  R)4x4 


(109s) 
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(109b) 


---R  ^  f 2(r;x)  a  1  +  (1  -  R)2x2  +  g  { (1  -  R)4  -  R[R2  +  (1  -  R)J}x**  . 


From  these  we  observe  that  the  tangent  is  increasing  faster  than  f2  and  hence 
will  continue  to  lie  above  f2  until  it  reaches  its  first  singular  point, 
ir/2(1  -  R).  It  will  then  go  negative  and  there  will  be  no  intersection  out  to 
x  *  x_(R).  Over  this  range  of  x,  f2(R;x)  is  positive  and  increasing  until  it 
goes  to  at  x_(R).  For  those  values  of  R  4  RQ,  such  that  x+(R)  <  3ir/2(1  - 


R).  there  will  be  no  intersection  for  x  £  x+(R).  The  first  intersection 
X21  »R>  will  occur  just  before  x  =  3n/2(1  -  R),  where  both  f2  and  tan(1  - 

D  W  Aro  nno i  ti  tto  .  Pnnaar«iif  i  in  fmaoi  n/to  »[  1  )  (Dl  will  fhan 


R)x  are  positive.  Consecutive  crossings  xii/(R)  will  then  occur  on  the  posi¬ 


tive-going  arcs  of  the  tangent  curve  with  almost  the  periodicity  of  the  tan¬ 
gent  function  itself.  (Fig.  17,  shown  later,  is  a  good  illustration  of  this 
behavior.)  For  larger  values  of  R  <  RQ,  when  x+(R)  occurs  for  x  >  3it/2(1  - 
R),  the  first  crossing  x^j)(R)  will  occur  on  the  negative  arc  of  the  first 
branch  of  the  tangent  after  its  second  singularity.  Successive  crossings  will 
occur  on  successive  negative  arcs  of  the  tangent  out  to  x+(R),  since  f2(R;x) 
is  negative  between  x_(R)  and  x+(R).  Thereafter  the  curves  will  intersect  on 
the  positive -going  branches  of  the  tangent  curve  with  nearly  the  periodicity 
of  the  tangent  function.  These  intersections  spread  apart  with  increasing  R, 
as  is  clearly  illustrated  in  figures  15  and  16. 


Figures  15  through  24  show  the  details  of  the  f2(R;x),  tan(1  -  R)x 
and  the  location  of  the  eigenvalues  *2p  (R)  for  x  $  20  for  R  -  0.05  and  for 
0.1  3  R  £  0.9  in  steps  of  0.1.  Although  these  illustrations  only  cover  cases 
(a)  and  (c),  the  above  discussions  serve  to  round  out  the  entire  picture  for 
determining  the  TE,  n  -  2,  eigenvalues. 
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Figure 

^(1) 

T 


1 7.  Generating  functions  for 
with  R  *  0.20. 


Figure  18.  Generating  functions  for 
yil*  with  R  =  0.30. 
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which  give,  upon  substitution  into  equation  (94),  and  some  simplification. 


0  = 


■J([^?)  ^  **(?**)  “®'?) 


cos  y 


♦  (4  - 1) 


V/ 


sm  y 


] 


.  [715  .  6_ V  sin  .  /I5  .  n  cos  ir/11  _  6  \  co3  x  +  /15  .  n  sin  1) 

L\y  y2)  [y  7  Jl>4  *2/  l*3  7  J] 


Again  we  see  that  the  real  part  of  the  expression  vanishes  identically.  As 
before,  we  perform  some  algebraic  rearranging  of  this  relation  to  the  form 


cos ( x  -  y ) 

or  again  to  the  preferred  form, 
tan[(1  -  R)x]  =  lx 


I1 5  -  r2*2Y1 

tJi  + 

P5  -  ^1 

p5  -  R2x2" 

|\l 

\5  -  2R2x2/_T 

V  9 

[5  -  2x2J 

L5  -  2R2x2. 

I/)- 

(111) 


values  X3p^(R)  at  which  intersection  of  the  functions  tan(1  -  R)x  and  f3(R;x) 
occur.  J 


Once  again  the  eigenvalue  problem  is  reduced  to  determining  the 


f 3(R;x) 


We  can  write  fj(R;x)  as  a  ratio  of  two  polynomials: 
3(1  -  R)x[  75  -  5(l  -  5R  +  R2)x2  +  2R2x**] 


225  -  45(2  -  5R  +  2R2JX2  -  3r(  5  -  12R  +  5R2)x4  +  R3x6 


(112) 


The  denominator  considered  as  a  polynomial  in  x2  indicates  that  the 
maximum  number  of  positive  roots  is  two  for  x2.  This  in  turn  tells  us  that 
there  can  be  at  most  two  positive  real  roots,  x_  and  x+.  For  small  values  of 
R  this  is  indeed  the  case.  At  these  values  x^,  f^  becomes  singular.  Since 
the  numerator  and  denominator  are  positive  for  small  x  and  very  large  x,  f3 
will  also  be  positive  for  these  ranges.  Just  as  in  the  TE  case  where  n  =  2, 


as  R  increases,  x+  and  x_  move  toward  one  another  and  at  some  value  R  =*  RQ  we 


have  a  single  real  positive  double  root  in  the  denominator.  ttius,  we  will 
have  a  considerable  amount  of  behavior  of  n  **  3  that  closely  resembles  that 
for  n  «  2.  Here  we  have  an  additional  characteristic,  namely  a  numerator  in 
equation  (112)  that  can  have  at  most  two  zeros  for  x  >  0.  This  will  somewhat 


39 


• 


m 


v. 


m 


m 


fc*1 

> 


£ 

V**  I 


& 


as 

S 


\v^ 
v7« 

‘.“A 


\\S| 


•  m\m  ■  «  */  ’/  . 


V 


modify  the  n  =  2  type  of  behavior.  In  figures  25  through  37  we  illustrate  the 
characteristics  of  f3(R;x),  tan(1  -  R)x,  and  the  intersections  xP '  (R)  <  20. 
Thus,  figure  25  shows  the  behavior  for  the  small  value  of  R  =  0.05.  We  can 
see  in  this  illustration  that  f3(R  =  0. 05; x)  has  a  zero  x_(0. 05)  in  the  denom¬ 
inator  below  x  =  2,  and  f3  is  very  close  to  tan(1  -  R)x  all  the  way  out  to 
about  ir/(  1  -  R). 


60.00 

Z  40.00 

R  = 

— i 

=  0 

1.05 

| _ 1 

i — 

UH 

| 

| 

O 

^  2000 

■ 

| 

£ 

1  0. 
1 

85  -20.00 
z 

UI 

°  -40  .00 

I 

1 

1 

□ 

fln 

E 

■ 

is 

n 

i  i 

i 

n 

m 

a  i 

_ i 

o - 

c 

3 

i 

i 

n 

■ 

m 

jl _ -I  i  I _ i  I  I _ I _ Li _ I 

0.  2.50  5.00  7.50  10.00  12.50  15.00  17.50  20.00 

X 


Figure  25.  Generating  functions  for 
Y  3p^  with  R  *  0.05. 


Figure  27.  Generating  functions  for 
Y3p^  with  R  =  0.10. 
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Figure  26.  Enlarged  view  of  generating 
functions  for  Y 3^  (0.05)  near  their 
first  singularities. 
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Figure  28.  Enlarc 
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trged  view  of  generating 
functions  for  yP'  (0.10)  near  their 
first  singularities. 


Figure  29.  Enlarged  view  of 
generating  functions  for 
y^p  (0.10)  showing  no 
intersections  just  beyond  first 
singularities. 
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for  R  =  0.10  since  the  tangent  has  been  stretched  out  along  x.  After  the 
second  singularity  in  f 3  is  passed,  the  intersections  occur  along  the  posi¬ 
tive-going  arcs  of  the  tangent  just  before  the  tangent  goes  singular. 

Increasing  R  still  further  causes  the  two  singular  points  of 
f3(R;x)  to  move  further  toward  each  other,  as  can  be  clearly  seen  in  figures 
31  and  32.  In  the  latter  case,  where  R  =  0.40,  we  note  that  the  two  roots  of 
the  denominator  of  f3(0.40)  lie  between  the  first  two  singularities  of  the 
tangent;  hence,  the  first  intersection  occurs  beyond  the  higher  value  of  x 
where  f3  goes  singular.  The  eigenvalues  behave  now  in  a  familiar  fashion. 

Further  increasing  R  eventually  causes  the  two  singularities  of  f3 
to  merge  into  one.  Thereafter,  increasing  R,  as  can  be  seen  in  figures  33 
through  37,  the  singularities  are  gone  from  f3(R;x)  and  we  again  have  the  same 
characteristics  that  we  found  for  f^R,^)  and  for  the  eigenvalues  in  TE,  for  n 
=  1. 

3.2.4  The  y£L*(R) 

The  last  TE  order  we  will  consider  explicitly  (and  thus  only  in  a 
limited  fashion)  is  that  for  n  =  4.  Here  we  have  for  the  pertinent  modified 
spherical  Bessel  functions. 


i4( jx) 


*  7^5-  *3=] 

(m .  .  i\  si„  x .  list .  i°\ . 

\*5  *3  v  \«*  *7 


— 1 
(jx)2J 


cosh( jx) 


(113) 


COS  X 


*4(jx)  =  e 


e”jx  f  105  105  45  10  _l1 

L(j*)»  (jx)4  (jx)3  (jx)2  3xJ 

1/105  10\  /l  05  45  .  1  \  .  1 

I-?*?  J 

K105  45  A  l\  /l05  10\  1 

""x3"  ~  x3"  +  COS  *  +  yx4"  ”  x2/  S^n  j  * 


(114) 


Substituting  into  equation  (94)  we  get 


0  -  j 


K1 05  45  l\  .  /l05  10\  ir/105  45  ,  l\ 

x5  *  x3  *)  8  n  *  "  (  x4  '  x2)  C°  J  [  \  y5  '  y3  y)  C°B  7 

/l  05  10\  .  1  r/105  45  l\  ,  /l  05  10\  1 

(lw)8inyJ  W‘?  +  i)8iny“  (Tw)C08yJ 

|7l05  45  l\  / 105  10\  .  1 

lU  -  J  +  x)  008  *  +  '  ?)  9in  *J  ' 


Once  aore  we  observe  that  the  real  part  vanishes  identically.  This  can  be 
rearranged  to 


and  this  in  turn  can  be  brought  into  the  form 
tan[(1  -  R)x] 


=  I  [7  105  ~  10x2  \  _  /  IPS  -  10R2x2  VI 

"  \%0S  -  45x2  +  x4/  \105  -  45R2x2  +  R4x4/J 

*  [l  +  Rx2 ( 1 05  -  1 0x2 ) ( 1 05  -  1 0R2 x2 )  +  ( 1 05  -  45x2  +  x4)(l05  -  45R2x2  +  R4x4)J 


=  f 4<R;x) 


(115) 


which  upon  simplification  is  a  rational  fraction  whose  numerator  is  a  polyno¬ 
mial  of  degree  7  whose  denominator  is  a  polynomial  of  fourth  degree  in  x2. 
The  coefficients  in  these  polynomials  are  of  course  functions  (i.e.,  polynomi¬ 
als)  of  the  ratio  of  inner  to  outer  spherical  radii  R. 

For  small  values  of  R,  the  denominator  will  have  four  real  positive 
roots.  If  we  examine  figure  39  for  R  =  0.05  we  can  see  the  two  smaller  of 
these  which  correspond  to  singularities  of  f^lRjx).  In  addition,  the  numera¬ 
tor  goes  to  zero  at  four  real  positive  values  of  x.  This  occurs  for  R  < 
0.05.  As  R  increases,  the  two  zeros  lying  between  the  inner  singularities  of 
f4(R;x)  merge,  and  by  R  =  0.05  they  are  gone,  leaving  only  the  innermost  and 
outermost  zeros  of  f^.  The  two  functions  f^(R;x)  and  tan(1  -  R)x  lie  ex¬ 
tremely  close  to  one  another  for  x  values  out  to  the  first  positive  zero  of  f  4 
as  can  be  seen  in  figures  38  through  43.  The  latter  of  these  is  for  R  = 
0.40.  Figure  39  shows  a  magnified  view  in  the  neighborhood  of  the  first 
singularities  for  both  functions  with  R  =>  0.05,  which  shows  rather  dramatic¬ 
ally  how  close  the  two  functions  are  indeed.  Fortunately,  this  figure  also 
shows  f4(0.05)  lying  below  tan  0. 95x  as  x  increases  from  zero,  the  tangent 
going  singular  first  and  the  relative  positions  of  the  two  functions  reversing 
as  x  increases  beyond  ttie  singularities.  Then,  in  figure  38,  we  can  note  that 
the  first  crossing  x^]'(0.05)  occurs  beyond  x  *  7.50  just  before  the  tangent 
goes  to  +*.  The  remaining  eigenvalues  for  x  £  20,  x^' (0.05),  occur  one  after 
the  other  on  the  negative-going  arcs  of  the  tangent,  as  can  be  very  plainly 
seen  in  the  figure. 
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As  R  increases,  the  outer  pairs  of  singularities  in  merge  into 
one  another,  which  has  already  taken  place  by  the  time  R  =  0.1.  Now  we  fur¬ 
ther  find  that  f4(0. 1;x)  again  has  two  zeros  at  positive  values  of  x.  This 
can  be  seen  in  figures  41  and  42  for  R  *  0.2  and  R  =  0.3. 


TE,  n 


Further  increase  in  R  merely  reproduces  the  behavior  we  found  for 
3.  This  can  be  readily  observed  by  examining  figures  43  through  47. 
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46.  Generating  functions  for  Figure  47.  Generating  functions  for 


with  R  ■  0.70. 


Y^p*  with  R  =  0.80. 


3.2.5  TE,  n  >  4 


Although  we  shall  not  go  into  any  detail  concerning  n  >  4,  we  can 
anticipate  the  additional  complexities  that  will  be  incurred  for  these  higher 
order  TE  mode  eigenvalues.  That  is,  additional  pairs  of  singularities  corre¬ 
sponding  to  a  larger  number  of  zeros  of  the  denominator  will  enter  and  de¬ 
crease  in  number  as  R  increases.  Similar  behavior  will  occur  in  the  fn(R;x) 
numerator.  Nevertheless,  fn  will  always  pass  through  the  origin  and  have  the 
same  slope  there  as  tan(1  -  R)x.  The  eigenvalues  will  always  be  such  that 
they  can  be  determined  in  the  manner  we  have  already  seen  from  the  roots  of 


tan(1  -  R)x  =  fn(R;x).  Even  though  the  formal  algebraic  analysis  for  the 
higher  n  can  be  quite  formidable,  it  can  be  carried  through  by  carefully  going 
back  and  forth  between  computer  calculations  and  graphics  and  formal  deriva¬ 
tions.  In  this  manner  all  singularities  and  zeros  of  the  fn(R;x)  can  be 
located,  as  well  as  the  sign  and  behavior  of  the  function  itself  over  perti¬ 
nent  ranges  of  R  and  x.  And,  of  course,  the  intersection  with  tan(1  -  R)x  can 
be  determined  definitively,  and  we  thereby  obtain  the  values.  In 

fact,  we  have  freely  used  this  approach  in  the  above  considerations  and  in 
those  below  for  the  TM  eigenvalues  Y^p^R)  wbich  we  shall  discuss  next. 

3.3  The  TM  Eigenvalues 

He  proceed  to  the  analysis  used  in  determining  the  TM  eigenvalues. 
The  roots  of  equation  (79)  have  to  be  found.  For  convenience,  we  repeat  the 
relation  here, 

[yain(Ya)] '[Ybkn(Yb] '  -  [Ybin(Yb)j '[Yakn(Ya)J '  =  0  , 


and  rewrite  it 


[ jxin(jx)] '[ jykn(jy)] '  -  [ jyin(jy)] '[ jxkn(jx)] •  =  0 


(116) 


As  in  the  TE  case,  due  to  the  nature  of  the  modified  spherical  Bessel  func¬ 
tions  there  will  be  roots  only  along  the  imaginary  axis,  and  they  will  occur 
as  complex  conjugate  pairs.  Only  those  along  the  positive  half  of  the  imagin¬ 
ary  axis  will  correspond  to  physically  acceptable  solutions  for  the  modes  of 
our  concentric  spherical  cavity. 

3. 3. 1  The  Yjp*(R) 

The  derivatives  that  we  shall  need  are 


[rakl<Y»)j;„jx  = 


[’  *  7*o*|  8inh 
j[(’  '  ?)  si" x 


sinh(jx)  - 


cosh( jx) 


(117) 


r 1  + 1  +  ii  - 

7l_  . 

\  ,  sin  x| 

L(jx)2  3*  J 

\x2 

/  x  J 

(118) 


..  •  .  ••  -  *  *  •  ■  • •  • ;  •  *  .  v .  * .  • 


.  .%  /  yv -v v v  v- *y 
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which  upon  substitution  in  equation  (116)  gives  us  after  some  rearranging 


-jk'-a-H'-M] 


sin(x  -  y) 


cos(x  -  y) 


(119) 


or,  equivalently, 

(1  -  R)x( 1  +  Rx2) 

tan((1  -  R)x]  *  — -  =  gj(R;x) 

[l  -(l  -R+R2)x2+  R2x4] 

Clearly  this  function,  gj(R;x),  is  of  the  same  form  as  f2(R;x)  that 
we  found  for  the  TE,  n  =  2  case.  We  have  two  zeros  in  the  denominator  and 
hence  two  distinct  singularities  for  g^R^)  for  R  <  RQ  *  (3  -  /5)/2  =■ 
0.381966.  These  singularities  occur  at 

/I  -  R  +  R2  .  /l  -  R  +  R2  r 

- *  /-^5 - P  • 


Also  in  this  case  g1  is  non-negative  for  0  £  x  <  x_;  it  is  negative  for 
x_  <  x  <  x+  and  positive  thereafter,  going  to  zero  asymptotically  as  x  becomes 
infinitely  large.  For  R  *  RQ,  the  two  singularities  merge  into  one  and 
gi(R0;x)  starts  at  zero  for  x  *  0  and  increases  monotonically  without  bound  as 
x  goes  to  the  double  root  at  x  =  ^3  +  /T/2  of  the  denominator.  On  the  other 
side  of  this  double  root,  g1(RQ;x)  decreases  monotonically  from  +»  to  zero  as 
x  increases.  Finally,  for  R  >  RQ  the  coefficient  of  x  in  the  denominator  of 
g.|(R;x)  is  positive  and  hence  g^(R;x)  is  well-behaved  starting  at  zero  at  the 
origin,  increasing  smoothly  to  a  maximum,  and  thereafter  decreasing  smoothly 
to  vanish  at  infinity,  remaining  non-negative  all  along  the  positive  side  of 
the  x-axis. 

For  R  <  RQ,  the  function  g^R.-x)  lies  above  the  tangent  function 
initially  and  goes  singular  first.  Then  it  is  negative  while  the  tangent  is 
still  positive  and  increasing.  Thus,  g^Rjx)  rises  from  earlier  than  the 
tangent,  after  x  =  ir/2{1  -  R),  and  still  lies  above  the  tangent  curve  for  a 
while.  Since  the  tangent  will  go  positive  at  x  *  ir/(1  -  R),  the  two  functions 
will  cross  before  tan(1  -  R)x  =  0.  This  behavior  will  persist  until  the 
second  singularity  of  g^(R;x)  occurs  before  x  =  w/2 ( 1  -  R).  We  see  this  in 
figures  48  through  51  for  R  =  0.05,  0.1,  0.15,  0.20,  and  0.30.  In  figure  52, 


vv> 


for  R  =  0.35,  the  second  singularity  occurs  before  ir/2(1  -  R)  and  hence  the 
first  intersection,  x{^(0.35),  now  occurs  while  the  tangent  is  going  to  its 
first  singularity  and  g^  is  returning  from  its  second  singularity.  In  figures 
48  through  51,  the  earlier  intersections  occur  on  the  negative  arcs  of  the 
tangent  curve  that  lie  between  the  singularities  of  g^(R;x).  Ihe  later  inter¬ 
sections  must  then  lie  on  the  positive  increasing  arcs  of  the  tangent,  as  can 
be  seen  in  these  same  figures.  Since  the  two  singularities  collapse  into  one 
at  R  *  R0  for  x  <  n/2(l  -  RQ) ,  we  can  see  in  figures  52  through  58  that  for 
r  £  R  intersections  or  eigenvalues  occur  on  each  successive  positive -going 
branch  of  the  tangent  curve.  As  R  increases,  the  successive  eigenvalues  move 
farther  apart  as  tan(1  -  R)x  spreads  out  along  x.  In  the  TM,  n  *  1  case,  we 
have  included  results  for  R  at  intermediate  values  to  our  more  or  less  stan¬ 
dard  interval  of  0.1  to  show  finer  detailed  behavior  of  the  indicator  function 
g1 (R;x). 
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Figure  48.  Generating  functions  for 
with  R  =  0.05. 
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with  R  =  0.90. 


3.3.2  The  yi^lR) 
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He  need  the  derivatives  of  the  products 
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.-3*1  6  +  _6  + 


$r«] 


\ 

■  [(?  - )  -  -  (?  - 1) si”  *]  -  {(?  -  ’) 8in  - + (?  -  j)  ^  *] 


(121 


Substituting  these  two  relations  into  equation  (116)  and  rearranging,  we  get 
0  *  j  {\l—  ~  — H -rr  ~  -1  +  (—  -  1  H  —  "  Ml  ain(x  -  y) 


or,  with  some  further  rearranging,  we  bring  this  to  the  desired  form: 


51 


tanl(1  -  R)x] 


3x[  ( 6  -  «?)(2  -  ##)  -  R(6  -  R2x2)(2  -  x2)] 
9(2  -  x^)(2  -  R^x^)  +  Rx2(6  -  R2x2) 


g2(R;x)  . 


(122) 


We  note  that  this  function  g2(R;x)  is  the  ratio  of  a  fourth -degree  polynomial 
to  a  sixth-degree  polynomial,  just  as  we  saw  earlier  for  the  TE,  n  =  3,  case. 
Actually,  the  denominator  of  g 2 ( R ; x )  is  a  cubic  in  x2.  Thus,  we  can  rewrite 
g2  as 


g2(R?x) 


3(1  -  R)x[  12  -  2(1  -  R)2x2  +  R2x4] 

36  -  18(1  -  R)2x2  -  3r[ 2(1  -  R)2  +  r]x4  +  R3x6 


(123) 


Examining  the  denominator  we  see  that  the  maximum  number  of  positive  roots  for 
x2  will  be  two,  and  in  turn  the  maximum  number  of  positive  roots  for  x  will  be 
two,  x_(R)  and  x+(R).  Now  let  us  look  at  the  numerator.  For  R  =  RQ  = 
(l  +  /3)  -  //3  +  2/~3  =  0.189,591,  there  will  be  only  one  root  for  the  numera¬ 
tor  at  xQ  =  x(rq)  =  (1  -  Rq)/R0  =  4.274,512.  For  R  <  RQ,  the  numerator  will 
have  two  roots  and  hence  g0  will  vanish  at  x  =  x„  , (R),  where 

**  Oft 

xo,±(R)  ±  /l  -  [  1 2R2 / ( 1  -  R)4  ]  j  , 


and 


xQ, _(R)  <  xQ  <  xo,  +  ^ ^  * 

Thus  we  see  that  as  R  increases,  xQ  _  and  xQ  +  move  together  and  coincide  at 
xQ  where  R  >  RQ.  For  R  >  RQ,  the  numerator  becomes  positive  definite. 

If  we  expand  g2(R;x)  around  x  =  0  in  a  power  series,  we  find  that, 
to  order  x5,  g2(R;x)  >  tan(1  -  R)x  near  the  origin.  Of  course  g2  and  the 
tangents  both  pass  through  the  origin  with  the  same  slope.  As  x  increases 
from  zero,  g2  will  have  a  larger  slope  than  the  tangent.  For  small  R,  the 
function  g2<R;x)  goes  singular  before  the  tangent  just  as  in  the  IM,  n  =  1, 
case.  As  a  consequence  of  this  behavior,  g2<R;x)  and  tan(1  -  R)x  will  not 
intersect  for  x  less  than  w/(1  -  R).  This  property  persists  until  R  is  suffi¬ 
ciently  large  for  g2(R;x)  to  have  no  roots  other  than  at  the  origin  or  infi¬ 
nity;  i.e.,  for  R  to  exceed  RQ.  These  characteristics  can  be  seen  in  figures 
59  through  62.  Although  explicit  analysis  can  be  carried  out  to  determine 
when  and  where,  as  a  function  of  R,  the  denominator  of  g2(R;x)  vanishes  and, 
hence,  g2  becomes  singular,  we  shall  not  do  so  here.  Examination  of  figures 
59  through  67  reveals  quite  readily  that,  as  R  increases,  there  are  two  singu¬ 
lar  values  of  x,  0  <  x_(R)  <  x+(R)  which  move  toward  each  other  as  R  increases 


52 


through  R  *  0. 58  and  which  coalesce  into  one  singularity  at  ti)e  double  root  of 
the  denominator.  By  R  =  0.60  this  singular  behavior  is  gone  and  as  we  see  in 
figures  67  through  71  we  have  a  familiar  behavior  for  g2(R»x).  For  those 
values  of  R  for  which  there  are  two  singularities,  g2(R?x)  is  a  monotonic 
increasing  function  greater  than  zero  and  increasing  with  x  increasing  to 
x_(R).  As  x  increases  further,  g2(R;x)  returns  from  toward  a  maximum  value 
(which  for  R  >  RQ  is  negative)  and  then  again  goes  monotonically  to  -*  as  x 
increases  toward  x+(R).  It  then  jumps  to  +»  and  decreases  monotonically  to  go 
asymptotically  to  zero  as  x  goes  to  infinity.  In  between  x_  and  x+  the  inter¬ 
sections  are  on  the  positive -going  arcs  of  the  tangent  for  g2  >  0  and  on  the 
negative  branches  for  g2  <  0.  For  x  >  x+  the  intersections  always  occur  on 
the  positive  arcs  of  the  tangent  function  and  are  very  close  to  being  periodic 
with  the  period  of  the  tangent  itself.  For  R  £  0.58,  the  first  crossing  of 


the  two  curves,  x^^(R),  occurs  on  the  first  positive  arc  of  the  tangent  (see 


fig.  67). 
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Figure  59.  Generating  functions  for  Figure  60.  Generating  functions  for 
with  R  =  0.05.  with  R  =  0.10. 
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61 .  Generating  functions  for 
with  R  =  0.15. 


Figure  62.  Generating  functions  for 
with  R  =  0.20. 
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GENERA 


[Yak3(Ya)]YaEjx 


_J1±_  +  -JLL-  +  _JLL_  +  ^  +  i 
L(jx)4  (jx)3  (jx)2  3X  J 

-(P  -%+  ')  “*  »  -  p  -  I)  8i"  »]  +  ip  -  I) 


(125) 


-6' ***’)]  coslx-,)}  '  ' 


Which  When  substituted  into  equation  (116)  gives,  after  some  rearranging. 


and  in  turn  this  can  be  put  in  the  desired  form  for  determining  the  eigen¬ 
values  from  the  intersections  of  two  curves,  namely. 


tan(1  -  R)x  =  3x 


15  -  2X2  _  /  15  -  2R2x2  \ 

[45  -  21  x2  +  x4  \45  -  21R2x2  +  R4x4/ 


(126) 


f  |l  +  9RX2 


(l5  -  2x2](l5  -  2R2x2 


(45  -  21X2  +  x4  )  (  45  -  21R2x2 


i _ 1) 

+  R4  x4  )J  j 


=  g3(R;x)  . 


We  point  out  that  the  denominator  of  g3(R;x)  is  a  quartic  in  x2, 
just  as  we  observed  for  the  polynomial  denominator  of  f4(R;x)  in  the  TE,  n  * 
4,  case.  Then  we  can  anticipate  a  rapid  increase  in  the  complexity  of  the 
extensive  analysis  that  accompanies  g3(R;x).  Only  a  small  portion  of  this 
analysis  is  included  here.  Thus,  for  the  behavior  near  x  *  0,  we  have 


g3(R?x)  5  3(1  -  R)x  [2R3x6  -  3r( 5R2  -  9R  +  5)x4  -  45( 2R2  -  5R  +  2)x2  +  675] 


*  [r4x®  -  3R2(  7R2  -  12R  +  7)x6  +  9(  5R4  -  30R3  +  49R2  -  30R  +  5)x4 

(127) 


-  1 45( 7R2  -  15R  +  7)x2  +  2025]  . 


Expanding  in  a  power  series  around  x  -  0,  we  find  to  order  x5 


g3(R;x) 


(1 


R)x 


(1 


R)3X3  2(1  -  R)5x5 

3  15 


9 


but  to  the  same  order. 


tan(1  -  R)x  "  (1  -  R)x  + 


(1  -  R) 3 x3 
3 


2(1  -  R)5x5 
15 


Thus  we  see  that  just  as  in  the  TE,  n  =  3,  case,  both  the  tangent  and  g3  coin¬ 
cide  to  order  x5.  This  is  reflected  in  the  two  curves  being  nearly  indistin¬ 
guishable  for  an  appreciable  range  in  x  for  the  smaller  values  of  R  of  figures 
72  through  80.  In  figure  72,  for  R  *  0.05,  in  fact,  the  two  functions  are 
extremely  close  together  even  beyond  the  first  singularity  of  g3(R;x)  and  the 
first  pole  of  tan(1  -  R)x.  Next,  if  we  examine  the  coefficient  in  the  denom¬ 
inator,  we  observe  that  for  small  R  there  are  at  most  four  real  positive  roots 
or  equivalently  at  most  four  singularities  in  g3(R;x)  for  positive  x.  The 
same  type  of  consideration  shows  that  the  numerator  has,  at  most,  two  roots 
for  small  R.  In  figure  72,  for  R  =  0.05,  we  see  one  of  the  zeros  of 
g3(0.05;x)  and  also  two  of  the  singularities.  For  x  beyond  the  range  shown  in 
this  figure,  the  pattern  is  repeated,  but  in  reverse  order;  i.e.,  as  x  in¬ 
creases  beyond  .20  we  find  g3(0.05;x)  negative  and  going  to  -<*>,  then  jumping  to 
•H»  and  decreasing  as  x  increases,  going  through  zero  at  a  specific  value  of  x 
and  then  going  to  again  as  x  increases  further.  Next,  g3  jumps  to  +■>  and 
then  goes  asymptotically  to  zero  as  x  As  R  increases,  the  singularities 
move  toward  each  other,  with  the  inner  pair  closing  together  more  rapidly. 
This  can  be  seen  in  figure  73  for  R  =  0.10.  Further  increasing  R  causes  the 
two  innermost  singularities  to  coincide  and  then  disappear  altogether.  When 
this  occurs,  g3(R;x)  has  two  zeros  between  the  remaining  singularities. 
Figures  74  and  75  show  this  rather  clearly.  Also,  we  note  that  the  outer 
singularities  move  inward  quite  rapidly  toward  the  slowly  changing  innermost 
singularity.  With  further  increase  in  R,  the  two  roots  of  g3(R;x)  move  to¬ 
gether,  coalesce,  and  then  disappear  altogether,  as  can  be  seen  in  figures  76, 
77,  and  78  for  R  *  0.40,  0.50,  and  0.60,  respectively.  Now  g2  is  less  than 
zero  in  between  the  singularities  and  positive  outside.  As  R  increases  even 
further,  the  two  remaining  singularities  move  together,  coalesce,  and  vanish. 
The  final  state  of  affairs  is  represented  by  figure  79  for  R  *  0. 70.  Further 
growth  in  the  size  of  R  as  usual  results  in  the  peak  of  g3(R;x)  dropping  and 
shifting  down  in  x.  All  the  features  of  g3(R;x)  found  here  for  TM,  n  «  3,  are 
now  familiar  properties. 
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functions  for 


Figure  81.  Generating  functions  for 
with  R  =  0.90. 


Next,  we  briefly  consider  the  intersections  (i.e.,  the  locations  of 
the  eigenvalues)  of  g3(R?x)  and  tan(1  -  R)x.  For  small  values  of  R  (see  fig. 
72  for  R  ■  0.05),  the  first  two  singularities  in  g3(R;x)  just  precede  those 
for  tan(1  -  R)x.  For  this  situation,  the  first  intersection  will  occur  just 
beyond  x  ■  3n/2(1  -  R),  i.e.,  on  the  negative  arcs  of  both  functions.  The 
next  few  eigenvalues  occur  between  the  inner  pair  of  singularities  on  the 
negative  arcs  of  the  tangent  just  after  the  tangent  returns  from  -°°.  The  bulk 
of  the  remaining  properties  of  the  locations  of  the  crossings  of  g3  with  the 
tangent  curve  has  already  been  described  in  the  foregoing  discussions.  These 
predicted  characteristics  can  be  verified  if  we  examine  the  corresponding  TM 
sets  of  curves  for  n  *  3  in  figures  72  through  81. 
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We  note  that  the  denominator  of  g4(R;x)  is  a  quintic  in  x2,  and  the 
numerator  is  the  product  of  x  and  a  quartic  in  x2.  This  is  a  level  of  com¬ 
plexity  beyond  those  met  thus  far,  but  the  general  nature  of  the  analysis 
required  is  pretty  much  the  same  as  that  we  have  carried  out  thus  far.  Upon 
examination  of  the  denominator  of  g^(R;x),  we  find  that  there  are  in  this  case 
at  most  four  positive  values  of  x  that  are  roots.  Hence,  g4(R;x)  will  have  at 
most  four  singularities  for  positive  x.  In  fact,  this  is  reflected  in  the 
behavior  of  g4(R;x)  as  shown  in  figures  82  through  84,  where  we  can  discern 
for  R  <  0.20  the  existence  of  four  singularities.  Figure  84  explicitly  shows 
for  R  =  0.20  three  of  the  four.  As  R  increases,  we  see  that  the  inner  pair 
move  together  to  coalesce  and  then  disappear  by  the  time  R  attains  the  value  R 
=  0.30  (see  fig.  85).  With  further  increase  in  R,  the  remaining  two  singular¬ 
ities  repeat  this  behavior,  and  by  R  =  0.80  all  the  singularity  in  g4(R;x) 
ceases.  For  R  >  0.80,  we  have  the  usual  well-behaved  smooth  increase  to  a 
maximum,  which  decreases  and  moves  inward  with  increasing  R.  This  behavior  of 
g4(R;x)  is  easily  seen  in  figures  90  and  91  for  R  =  0.80  and  0.90,  respec¬ 
tively.  The  numerator  of  g4(R;x)  is  a  product  of  x  with  a  polynomial  that  is 
a  quartic  in  x2.  Examination  of  this  numerator  shows  that  g4(R;x)  will 
possess  a  maximum  of  four  roots  for  positive  x.  For  small  R  these  are 
located,  as  we  have  seen  in  other  cases,  as  follows:  one  each  between  the 
inner  and  the  outer  pairs  of  singularities  and  two  between  the  inner  pair  of 
singularities.  This  can  be  seen  in  figure  82  for  R  =  0.05.  As  R  increases, 
the  innermost  pair  of  roots  move  together,  then  merge,  and  by  R  =  0.10  (see 
fig.  83),  disappear,  leaving  g4  negative  between  the  pair  of  inner  singular¬ 
ities.  As  R  increases  further,  the  two  inner  singularities  disappear,  and  two 
roots  are  present  between  the  remaining  pair  of  singularities  (see  fig.  85). 
As  R  grows  larger,  these  two  zeros  move  towards  one  another,  merge,  and  disap¬ 
pear.  For  R  =  0.60,  as  can  be  plainly  seen  in  figure  88,  there  are  no  finite 
positive  values  of  x  for  which  g4(R;x)  vanishes. 
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with  R  =  0.80. 


y£2)  with  R  =  0.90. 


For  future  reference,  we  write  out  g4(R;x)  in  reduced  form: 
g4(R;x)  =  (1  -  R)x[l  -  -g^n  -  28R  +  11R2)x2  +  357^0  (28  "  252R  +  463r2  '  25‘ 

+  28R4)x4  -  5R2(  39  -  71R  +  39R2)x6  +  10R4x8J  *  -  ( 1 3  -  28R  +  13R2)x2 


(131) 


+  (56  -  308R  +  507R2  -  30 8R3  +  56R4 )  x4 

+  -  \  —  r(  84  -  390R  +  605R2  -  390R3  +  84R4lx6 
35,280  t  J 

-  Is^So  ("  -  20R  ♦  ,,r2)xS  *-ni7S5o]  • 


As  can  be  seen  in  figures  82  through  85,  the  two  curves  lie  very  close  togeth¬ 
er  for  an  appreciable  range  of  x  when  R  <  0.30.  Careful  analysis  of  these 
expansions  shows  that  g4(R;x)  lies  above  the  tan(1  -  R)x.  For  the  small  R 
values,  g4(R;x)  attains  its  first  singularity  before  the  tangent  function. 
For  larger  R  this  characteristic  is  plainly  seen  (fig.  87  through  89  for  R  * 
0.50  to  0.70).  From  the  behavior  of  the  two  functions  for  the  smaller  values 
of  R  which  correspond  to  the  presence  of  singularities  in  g4(R;x),  it  is 
relatively  easy  to  discern  the  fact  that  there  will  be  no  intersection  until 
at  least  x  =  6.  However,  as  R  increases  beyond  r  =  0.40  (see  fig.  86,  for 
example),  the  first  eigenvalue  x^^(R)  occurs  at  a  smaller  value  of  x.  All 
the  remaining  features  of  the  behavior  of  the  intersection  positions  have 
essentially  been  discussed  in  the  foregoing  cases. 


3.3.5  The  Y*p  (R) 

For  the  final  case,  we  shall  present  an  explicit  discussion  for  TM 
5.  Our  two  starting  relations  are 


[Y«i, (T.)]'  *  fiHS  .  ZS!L  .  -U»  .  ,1  ,lnMjlt) 

ra-]x  L(  jx)6  (jx)4  (jx)2  J 

. riZ2§_ . -«sl .  »-\  cosh(jx)  $  Jfi ip. 


1315  -  42x2  + 
- - - 

4725  -  2205x2  +  i: 
R(  315  -  42R2x2  +  R4x4)  "1 

-  2205R2  x2  +  1 2 OR4  x4  -  R6x6J 


2  Ox4  -  x6 


(134) 


*  [l  +  225Rx2(  315  -  42X2  +  15x4)(315  -  42R2x2  +  15R4x4) 


*  (4725  -  2205X2  +  120x4  -  x6)(4725  -  2205R2x2  +  120R4x4  -  R6x6)]  =  g5(R;x) 


As  we  can  anticipate  at  this  stage  of  the  discussion,  g5(R;x)  is  a 
ratio  of  two  polynomials.  Its  denominator  is  a  sixth -degree  polynomial  in 
x2.  It  has  as  its  numerator  the  product  of  x  and  a  quintic  in  x2.  Again  it 
would  appear  that  we  have  to  deal  with  a  higher  degree  of  complexity  than 
experienced  earlier  in  the  properties  of  the  two  functions  that  generate  the 
eigenvalues.  For  reference  purposes  we  give  the  reduced  form  of  g5(R;x): 


g5(R;x)  S  (1  -  R)x[  1 , 460, 025  -  4,725(42  -  5R  +  42R2)x2  +  945(5  -  35R  +  63R2 
-  3SR3  +  5R4 }x4  +  1  5r(  2 j  -  1 24R  +  212R2  -  124R3  +  23R4)x6  -  6R3(7  -  13R  +  7R2)x8 


+  R4x10]  *  [1,460,025  -  99,225(7  -  15R  +  7R2)x2  +  945(40  -  210R  +  343R2  -  210R3 


+  40R4 ) x4  -  305(1  -  225R  +  56R2  -  84R3  +  56R4  -  225R5  +  R6)x6  +  3R2(49  -  3, 1 50R 


+  320R2  -  3,1 50R3  +  49R4)x8  -  R4(8  -  3,375R  +  8R2)x10  +  R6x12]  . 


Again  examining  the  numerator  of  g5(R;x),  we  observe  that  the 
maximum  number  of  finite  positive  roots  is  four.  As  usual,  g5  and  tan(1  -  R)x 
and  their  slopes  coincide  at  the  origin.  This  time,  the  denominator  shows  by 
simple  examination  that  the  maximum  number  of  its  positive  zeros  in  x2  is 
six.  However,  if  we  look  closely  at  the  coefficient  of  x10  in  the  denomina¬ 
tor,  we  observe  that  its  sign  will  be  positive  if  R  >  0.002370.  This  corre¬ 
sponds  to  an  extremely  small  inner  sphere.  For  R  values  exceeding  this  ratio 
just  cited,  there  will  be  at  most  four  singularities  for  g^CRjx)  on  the  posi¬ 
tive  x-axis.  At  this  time,  we  choose  to  ignore  the  characteristics  of  the 
extremely  small  R  case. 


Figures  92  through  95,  for  R  from  0.05  to  0.30,  display  the  situ¬ 
ation  for  the  range  of  R  where  ggfRf-x)  possesses  four  singularities  and  the 
number  of  zeros  goes  from  four  to  two.  Figures  96  through  99  (for  R  =  0.40  to 
R  =  0.70)  show  the  characteristics  of  g^(R;x)  for  the  range  of  R  where  there 
are  only  two  singularities  and  the  number  of  positive  finite  roots  goes  from 
two  to  zero. 


The  essential  behavior  of  g^(R;x),  tan(1  -  R)x,  and  the  eigenvalues 
is  very  much  like  that  we  have  already  seen  (except  for  R  <  0.002374).  Fig¬ 
ures  100  and  101  for  R  =  0.80  and  0.90,  respectively,  are  now  quite  familiar 
to  us  for  the  situation  where  g^(R;x)  has  no  singularities. 


3. 4  Summary  Discussion  of  Eigenvalues 


In  the  preceding  section  we  presented  in  various  degrees  of  extensive 
detail  the  behavior  of  the  appropriate  functions  whose  intersections  yield  the 
eigenvalues  for  a  large  number  of  the  lower-lying  TE  and  TM  modes.  We  observe 
that  as  the  order  index  n  increased,  invariably  the  functions  fn(R;x)  and 
gn(R;x)  became  more  complex.  Where  we  could,  we  merely  referred  to  common 
characteristics  that  appeared  for  lower  orders,  which  may  seem  to  imply  seri¬ 
ous  omission  of  detail  in  the  discussion.  Fortunately,  the  resulting  discus¬ 
sion  should  be  quite  adequate  when  one  properly  carries  over  the  portions  of 
earlier  detailed  discussions.  These  discussions  of  each  individual  case  and 
the  large  number  of  accompanying  illustrations  should  be  ample  enough  re¬ 
sources  to  establish  that  the  eigenvalues  can  be  determined  analytically  with 
an  investment  of  effort  of  modest  difficulty  at  best.  It  should  be  noted, 
however,  that  although  the  functions  fn  and  gn  may  present  themselves  as 
somewhat  formidable  obstacles  for  determining  the  eigenvalues,  this  is  cer¬ 
tainly  not  the  actual  situation,  except  for  a  small  number  of  eigenvalues. 
Even  this  small  set  can  be  made  tractable  with  additional  analytic  effort.  In 
general,  once  any  eigenvalue,  say  the  p^1,  Y^p  »  is  found  for  fixed  j  and  n, 
good  approximate  values  of  the  others  for  that  j,n  can  be  readily  determined. 
This  is  a  consequence  of  the  very  close  approximation  to  periodicity  that  the 
set  displays.  This  property  is  evident  in  figures  102  and  103.  Figure  102 
displays  the  manner  in  which  the  TM  eigenvalues  ( Y^pa)  vary  with  R,  the  ratio 
of  inner  to  outer  radii  for  the  non-lossy  dielectric,  where  n  goes  from  1  to 
5.  The  corresponding  situation  for  the  TE  eigenvalues  (Y^l^a)  for  n  from  1  to 
4  is  shown  in  figure  103.  In  the  interest  of  clarity,  these  sets  of  eigen¬ 
values  were  divided  up  and  displayed  as  shown  in  these  two  figures.  Most  of 
the  numerical  values  from  which  the  curves  were  obtained  are  given  in  table 
1 .  In  appendix  B  we  present  the  computer  program  that  generated  the  curves  in 
the  figures  and  the  values  in  table  1. 


We  have  already  mentioned  the  approximate  periodicity  with  p  for 
given  j  and  n  that  the  eigenvalues  display.  Next  we  point  out  that  the  proper 
behavior  is  displayed  by  the  Ypi  (R)  as  R  decreases;  i.e.,  they  go  over 
smoothly  to  the  R  =  0  or  empty  spherical  cavity  modes  given,  for  example,  by 
Harrington.4 


4J?.  F.  Harrington,  Time-Harmonic  Electromagnetic  Fields,  Tables  6,1  and  6,2, 
McGraw-Hill  Book  Co.,  Inc.,  NY  (1961),  270. 
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TABLE  1.  NUMBUCAL  VALUES  OP  THE  CONCENTRIC  SPHERICAL  CAVITY 
EIGENVALUES  y££*UU«  =  X*J'(R)  (FOR  VALUES  NOT  EXCEEDING  20) 


TABLE 

la.  THE  TS 

EIGENVALUES  x<’ * 
np 

(R) 

n 

m 

1 

P  -  1 

P  -  2 

P  -  3 

p  -  4 

p  -  5 

R 

m 

0.10 

4.52229 

7. 84656 

11.18355 

14.55528 

1 7. 95607 

R 

m 

0.20 

4.68640 

8.37806 

12.16588 

16.00839 

19.88013 

R 

m 

0. 30 

5.04272 

9.31415 

13.70082 

18.13297 

R 

m 

0.40 

5.63899 

10.69922 

15.86357 

R 

u 

0.50 

6.57200 

12.72134 

18. 95438 

R 

m 

0.60 

8.05525 

IS.81259 

R 

m 

0. 70 

10.60490 

R 

m 

0.80 

1 5. 78670 

R  -  0.05 

R  -  0.10 

R  -  0.20 

R  -  0.30 

R  -  0.40 

R  -  0.50 

R  -  0.60 

R  -  0.70 

R  -  0.80 

5.76350 

5. 76489 
5.79965 
5.96125 
6.35744 
7.11156 
8.44278 
10.86580 
15.94308 

9.09544 

9.10671 

9.31046 

9. 95545 
11.14065 
13.02612 
16.01949 

12.32478 

12.36709 

12.90198 

14. 16502 

16. 17099 
19.16249 

15.52005 
15.62634 
16.59799 
18. 49099 

n  -  3 

P  *  1 

p  -  2 

p  -  3 

p  *  4 

R  -  0.05 

6.98792 

10.41712 

13.69806 

16.92349 

R  -  0.10 

6. 98800 

10.41787 

13.70239 

16.93949 

R  -  0.20 

6.99344 

10.47287 

13.91692 

17.45049 

R  -  0.30 

7.05135 

10.83867 

14.83834 

19.01999 

R  -  0.40 

7.28038 

11.77298 

16.62299 

R  »  0.50 

7.84504 

13.47112 

19.47099 

R  •  0.60 

8.991 27 

16.32599 

R  »  0.70 

11.24SS6 

R  -  0.80 

16.17449 

18.70149 

18.90949 


n  ■  4 

P  »  1 

p  -  2 

p  -  3 

p  -  4 

R  -  0.05 

8.18256 

11.70491 

15.03967 

18. 30099 

R  «  0.10 

8.18256 

11.70495 

15.03998 

18. 30249 

R  -  0.20 

8.18328 

11.71655 

15.10656 

18.51699 

R  -  0.30 

8. 20089 

11.89008 

15.69443 

1 9. 70949 

R  -  0.40 

8.31928 

12.56665 

17.21049 

R  -  0.50 

8.71681 

14.04366 

19.  87549 

R  -  0.60 

9.67169 

16.72599 

R  -  0.70 

11.73234 

R  -  0.80 

16.47849 

/V 


R  -  0.05 

2.74162 

6.09682 

9.24655 

12.31915 

15.33189 

18.31399 

R  -  0.10 

2.72705 

5.96203 

8. 86383 

1 1 . 77078 

14.87325 

18.13634 

R  •  0.15 

2. 68855 

5.69396 

8.55252 

11.80250 

IS.  28791 

18. 86890 

R  -  0.20 

2.61986 

5.44453 

8.60724 

12.24786 

16.04S27 

1 9. 89957 

R  «  0.30 

2.41244 

5.34975 

9.  37465 

13.72043 

18.14151 

R  -  0.35 

2.29633 

5.49884 

9.97643 

14.70131 

19.4826S 

R  •  0.40 

2.18403 

5.75930 

10.71862 

15.86959 

R  «  0.50 

1 . 9846 

6.6185 

12.7281 

18.9565 

R  -  0.60 

1.82158 

8.07212 

15.81490 

1.68939 
1.58100 
1 . 49067 
1.42100 


10.61018 
1 5. 78800 


R  -  0.05 

4.97340 

8.72175 

12.06357 

15.31344 

18.52349 

R  -  0.10 

4. 97341 

8. 72142 

1 2. 06080 

15. 30030 

18.47949 

R  -  0.15 

4.97327 

8.71675 

12.02450 

15.15368 

18.12549 

R  -  0.20 

4. 97235 

8.68951 

11.85647 

14. 77432 

17.86149 

R  -  0.30 

4.95761 

8.41775 

11.37826 

15.00803 

19.08649 

R  -  0.40 

4.  881 38 

8.03405 

11.92239 

16.66549 

R  -  0.50 

4.68865 

8.16066 

13.51647 

19.48399 

R  «  0.60 

4.41209 

9.10303 

16. 34049 

R  -  0.70 

4.12700 

11.27912 

R  -  0.80 

3.87106 

16.18199 

R  -  0.90 

3.65137 

R  -  0.99 

3.  48150 
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TABLE  1.  NUMERICAL  VALUES  OF  THE  CONCENTRIC  SPHERICAL  CAVITY  EIGENVALUES 
=  Jt'i’lR)  (FOR  VALUES  NOT  EXCEEDING  20)  (Cont'd) 


np 

TABLE  1b. 


THE  TM  EIGENVALUES  (R) 
np 


6.06194 
6.06194 
6.06185 
6.05884 
6.03032 
5. 90808 
5.64693 
5.31651 
4. 99586 
4.71383 
4. 49450 


9. 96754 
9.96753 
9.96279 
9.85890 
9.46274 
9.25947 
9. 86902 
1 1 . 79056 
16.49099 


13.38012 
13.37997 
13.33330 
12. 86632 
12.86093 
14.12652 
16.75099 


16.67399 
16.67299 
16.46049 
16.05899 
1 7. 29099 
19.89849 


19. 91499 
19.91099 
1 9. 40899 
1 9. 84549 


The  variation  of  the  eigenvalues  with  R  displays  a  very  rich  diver¬ 
sity  of  characteristics.  We  shall  emphasize  only  a  restricted  number  of 
these,  which  is  far  from  exhaustive.  One  immediately  evident  characteristic 
is  the  grouping  of  the  eigenvalues  into  families  of  the  form  Y^p  (R)  for  j  and 
p  fixed  and  the  index  n  varying.  If  we  fix  our  attention  on  the  spectral 
group  Y^]  (R)a,  we  observe  for  these  TE  eigenvalues  that  as  R  increases  from 
zero  the  eigenvalues  vary  quite  slowly  at  first.  The  larger  the  value  of  the 
index  n,  the  greater  the  range  of  R  with  this  slow  variation.  This  behavior 
is  repeated  for  the  family  belonging  to  p  *  2.  Now,  however,  there  is  a 
smaller  range  of  R  over  which  the  variation  in  Y^2  lR)a  is  slow.  As  p  in¬ 
creases,  we  observe  in  figure  103  the  continuation  of  this  trend.  Immediately 
following  this  very  slow  variation  of  eigenvalue  with  increasing  R,  we  find  in 
each  case  in  figure  103  a  region  we  call  the  "knee"  over  which  this  rate  of 
change  speeds  up  from  slow  to  rapid  variation  with  R.  Once  past  the  knee  of 
the  curve,  Y^l^a  increases  much  more  slowly  as  R  increases.  What  we  have  here 
is  a  rough  qualitative  division  of  the  concentric  spherical  cavity  TE  eigen¬ 
values  into  three  regimes  corresponding  to  a  small  interior  sphere,  an  inter¬ 
mediate  sized  interior  sphere,  and  lastly  a  large  enclosed  sphere.  Before 
looking  at  figure  102  for  the  TM  eigenvalues,  we  note  that  all  the  Y„p  (R) 
always  increase  with  R  whether  variation  is  rapid  or  slow. 

Now  let  us  examine  the  manner  in  which  the  TM  eigenvalues  vary  as  R 
is  varied.  The  same  general  characteristics  are  displayed  by  the  Y^p*  for 
p  ^  2,  i.e.,  they  form  distinct  families  for  each  p?  also  they  indicate  that 


i 


three  eigenvalue  regimes  can  be  identified  which  correspond  to  small,  inter¬ 
mediate,  or  large  inner  spheres.  Thus  we  see,  for  the  p  >  2  families,  at 
first  a  slow  change  in  yj,? a  as  R  increases  from  zero.  The  larger  n  is,  the 
larger  is  this  range  of  R.  This  is  followed  by  the  intermediate  region  in 
which  this  variation  goes  from  slow  to  quite  rapid  and  then  in  turn  the  final 
region  of  rapid  variation  in  with  R.  In  the  TM  case  we  have  a  different 

type  of  behavior  that  is  quite  evident  in  figure  102.  Unlike  the  simple  knee 
of  the  TE  eigenvalue  trajectory,  we  now  have  a  "nose"  in  the  intermediate  TM 
region.  This  feature  indicates  the  presence  of  an  interesting  property, 
namely  that  for  p  >  2  as  the  inner  sphere  radius  grows  in  size  with  the  outer 

radius  fixed,  the  TM  eigenfrequency  decreases  initially  (for  a  not  insignifi¬ 

cant  range  of  R  in  many  cases )  and  then  later  increases  with  R  as  one  would 
expect.  Hie  behavior  for  large  R  is  that  where  the  inner  sphere  is  quite  close 

to  the  outer  sphere.  It  would  seem  to  be  quite  natural  in  this  situation  that 

the  eigenf requencies  be  such  that  multiples  of  half-wavelengths  span  the  gap 
between  inner  and  outer  spheres.  For  the  TE  modes  (see  fig.  103)  and  for  the 
p  >  2  TM  modes  (see  fig.  102)  this  is  indeed  the  case. 

For  R  large,  the  conducting  boundary  surfaces  are  close  together  and, 
parallel.  In  any  small  region  containing  portions  of  these  boundary  surfaces, 
we  have  what  approximately  looks  like  a  parallel-plate  capacitor  system.  Me 
should  have  some  natural  field  distribution  that  looks  like  the  parallel-plate 
system.  The  TE  solutions  cannot  satisfy  this.  What  we  do  observe  is  that  the 
TM  family  of  Y^?^R)  values  fills  the  bill  here.  Each  trajectory  in  this 
family  starts  out  at  the  empty  spherical  cavity  eigenvalue  (i.e.,  for  R  *  0) 
and  then,  as  R  increases,  decreases  in  the  familiar  manner  initially.  Now 
however,  continues  to  decrease  as  R  continues  to  increase.  This 

behavior  is  strikingly  unlike  all  the  other  trajectories  of  eigenvalues. 

Another  feature  of  considerable  significance  is  readily  evident  in 
figure  102  for  the  TM  eigenvalues  and  also  in  figure  103  for  the  TE  eigen¬ 
values.  Whenever  the  trajectories  intersect  we  have  degeneracy  present.  This 
degeneracy  is  present  in  addition  to  the  usual  degeneracy  associated  with  the 
angular  variables  0  and  $.  It  should  be  recalled  that  this  property  does  not 
exist  for  the  empty  spherical  cavity.  For  the  systems  that  we  shall  be  dis¬ 
cussing  in  subsequent  reports,  that  are  derived  from  the  concentric  spherical 
cavity  problem,  this  feature  shall  play  a  very  important  role.  Note  that  for 
a  given  value  of  R  we  observe  at  most  a  double  degeneracy  occurring. 

Some  further  interesting  features  that  can  be  extracted  from  figures  102 
and  103  are  the  ordering  of  the  eigenvalues.  Clearly,  for  very  small  R 
values,  the  ordering  is  precisely  that  of  the  empty  spherical  cavity.  As  R 
increases,  we  can  observe  switching  of  order  of  pairs  of  eigenvalues  as  the 
trajectories  cross.  In  addition  to  this  rearranging  of  the  ordering  of  the 
eigenvalues  due  to  degeneracy,  another  contributing  factor  enters.  This  is 
the  result  of  the  "nose"  region  in  the  TM  trajectories  which  move  the  eigen¬ 
values  to  lower  values  with  increasing  R;  e.g.,  at  R  •  0.3  the  TM  value  of 
lies  lower  than  the  TE  value  of  y^]'(0.3). 
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Although  we  have  explici tly  evaluated  a  rather  large  number  of  eigen¬ 
values  yII  (R)a  up  to  the  value  20,  we  have  by  no  means  determined  them  all. 


He  do  have,  however,  all  those  for  values  up  to  about  6.50. 
density  of  eigenvalues  increases  rather  quickly. 


Clearly,  the 


Further  examination  of  the  eigenvalue  trajectories  will  reveal  many 
other  physically  important  features.  Since  each  of  these  will  be  dealt  with 
in  subsequent  reports,  we  shall  defer  any  further  discussion  until  they 
appear. 


The  cavity  field  distributions  have  been  obtained  for  each  of  the 
modes  in  table  1.  Discussion  in  depth  of  the  field  details  constitutes  the 
next  report  in  the  sequence. 


APPENDIX  A 


Combining  the  two  pieces  we  obtain  after  integrating  over  $ 


n6 


,  /^  dr  r2  ♦np(Yip)r)^n*pt(Yi'p'r)  /*  <*9  sin  0  .  P^lcos  e)Pn.m( 


mm 


(cos  9) 


dPnm(cos  0)  dPn.m(cos  0) 

+  - - - -  • 

d0  d0 

However,  the  integral  over  0  is  given  in  Stratton1  as 


2n(n  +  1)  (n  +  m)l 


Jnn ' 


2n  +  1  (n  -  m)i 
We  then  have  remaining  integration  over  the  radius,  namely. 


2n(n  +  1)  (n  -  m)l  r  ^ 

w^nn' ^mm'  2n  +  "i —  (n  +  m)  t  d  ^  ♦npUnp  rJ4>np.  lYnp'r ) 

Now  4>  satisfies  the  differential  equation 

T  nD  v  J  *  nD  ' 


(i) 


rnp^ J  1  np 


o  -  (jYnp)r)2  — - ♦np(jYnp,r)  +  ^j^ip^)  (1)  ♦nptiYnp  r) 

d(jYAp *)2  d(jYnp  r) 

-  [(jYnp>r)2  +  n<n  +  1  ^♦np(^ip)r)  , 


and  ♦np,(3Ynp'r)  satisfies  the  equation 

0  -  (j^ip-r)2  (f}  2-»np’(Jirnp'r)  +  2(jyiplr) 

d(jYnp'rr 

-  [(j^ip'r)2  +  n(n  +  \  )]♦„?'(  jKnp*r)  • 


Then  multiplying  the  first  of  these  by  $np»(jYnp'r)  and  ^  second  by 

4  flY^'rl  and  subtracting  will  give  us,  after  some  further  reduction, 

T  npv  J  np  ' 


1 j.  Stratton,  Electromagnetic  Theory,  McGraw-Hill  Book  Co.,  Inc.,  NY  (1941), 
417. 


♦npliYnp  f)  _12E_li_2E  jj  ,  _ 


d$np( 3Ynp  r) 


(1)  ,  d$np' ( jYnp'r) 


The  first  term  on  the  right-hand  side  vanishes  since  the  ♦np(jT^p  r)  vanish  at 
the  cavity  boundaries.  The  remaining  term  cancels  the  second  pair  of  terms  in 


our  original  integral  over  r.  Thus,  we  have  left 


.H)^2  rb  _2  x  rx„n) 


[(r^')2  -  (yApO2]  I"  r2  -  0  . 

£1 

For  p'  *  p  the  scalar  functions  are  orthogonal  and  hence 

ssi  *(?)*£;0,(J)  •  -  [<;0>]2«nn'«pp'«»-  • 

cone 

spher 

cavity 

Next  we  consider  the  TE  even  modes.  By  following  the  same  steps  we  find 
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|“7  [  r^np(  3Ynp>r)  ]  -  [Y^>]2[r*np(jY^,r)]|  {-£  [r*n.p,(  jY^.r)] 

-  [Yn^p']2[r4'n1p'(jYn^p'0]J 

+  {r  I?  [r*np(^ip)r)])  (7  |f  [^n'p'CjY^p.r)]! 


pnm^cos  ^^n'm'^008  ®)  sin  sin 

dPnm^cos  0) 


dPntm* (cos  8 ) 
d9 


sin  m$  sin  m'4  +  mm1 


7 17  [r*np(jYip)r)j] 


7  5F  tr*n'p'(iYn'p*r)]J 


*  pnm^cos  0)  ^'m'^008  ®)  cos  cos 


Integration  over  <p  gives  tt6  ,  for  each  of  the  first  two  terms  and  _ ,  for 

mui  nun 

the  third  term.  In  turn,  integration  over  8  of  the  sum  of  the  second  and 
third  terms  gives 


2irn(n  +  1 )  (n  +  m)t  , 

2n  +  1  (n  -  m ) 1  ^nn'^am’  ' 

while  integration  over  8  after  $  gives 

2n  (n  +  m)  1 
2n  +  1  (n  -  m)I 


He  have  at  this  stage 


(jOJE)2  ///  dx 


>(2,o) 

“npm 


E 


(2,o) 

n'p'm' 


2ff  (n  +  m)  1  . 

2n  +  1  (n  -  m)l  °nn*»mm' 


X  ^ a  dr  .  {^T  IrWjYi>P>r)]  -  [Yip^]2  [r*np(^ip)r)]} 

*  {^T  tr*np*(jlTnp'r)]  -  [Yip-]2[r^np.(jriplr)]} 

+  n(n  +  1)  {7  aF  [r*np(^ip)r)]|  *  {?  d7  tr*np'(  jYnp?>r)]J  . 


If  we  use  the  differential  equation  that  defines  the  ^np( »  namely 
that  of  equation  (48),  we  can  rewrite  the  integral  on  the  rignt-hand  side  as 
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n(n  +  1)  /a  dr  {n(n  +  1  >*np(  jY^Mfnp'C jY^r)  +  |^  [rif>np( jr^r)] 

*  |j-  t^np'fiYnp'r)])  • 

Consider  the  integral  of  the  second  term  which  we  can  evaluate  as  follows: 


rb  .  d 


r  dr  §7  ['♦npCttiJM]  ^  [ripnp.(  jTTnp.r)]  = 

{^npCjYnp^)  ~  [r^np.(jyiplr)]J  ^ 

-  dr  ri|>np(  jYnp>r)  [ r^np.(  jTrip-r) ]  =  -  /^  dr  r<|>np(  jY^r) 

*  ([Yip-]2[ri|;np.(jYip]r)]  +  °(n^-1-)  [r»np.( jYnpM]}  . 

Combining  the  pieces  then,  we  have 

i  -imp  rrr  JT  i.(2,o)('  +  t  .  -2tn(n  +  1)(n  +■  m)l  f  .  r_  ( 2 )  -i  2 

(joje)  JJJ  dx  Enpm  L r J  En«pimi  *  Onn'^mm'LYnp'J 

(2n  +  l ) (n  -  in)] 

•  J**  dr  r2  •  <l'np(  i Ynp  ^ r ) ♦np ' ( jYnp ' r )  • 

Processing  this  remaining  integral  exactly  as  we  did  the  corresponding  one  for 

to  to  .  we  will  find  the  integral  over  the  radius  r  from  a  to  b  vanishes  for 
TnpTnp 

p  *  p ' .  Thus ,  we  have 


III  at 


PP'  * 


Next  we  demonstrate  orthogonality  between  the  Enpm°*  and  the  Enpm6 ^ ' 


Thus,  consider 


///  dx(r)!^e)(r)  •  *  1“  dr  r2^np(jTnp,r)^n'p‘(jYn1'p'r)  /’  d0  sin  6 


rb  _2>.  ra«(D_w  .  f6 


•  J^dtp  |  — .  -  •  mm ' Pnj, ( cos  0)Pn.m.(cos  0)  sin  m<f»  cos  m'$ 
a  [sinz0 

dPnatcos  0)  dPn.„. (cos  0)  \ 

+  — - - -  cos  »*  sin  «'♦)  . 
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For  precisely  the  same  reason  we  obtain 


III  *(?)  *&*’(*)  •  -  o 


Now  consider 


fff  t('-o)  .  t(2,e)  _  _J_  fb  2  f*  1  1  3lfrnpm 

JJJ  dT  Enpm  En'p’m1  “  “  jMe  dr  r  JQ  d*  (sin  0  3* 

3  / 1  3  r_,(e)  i\  1  8<^npm  3  /  1  3  r_,(e)  i\ 

30  (r  30  Lr*n'p'm'J)  *  sin  0  30  *  34,  \r  3r  Ir*n'p*m*  jj  = 


^-/b 

Ja 


/  dr  r2 


!q  d8  IQ  d4>  |«>np(^p)r)  7  I?  [rfn'p'Clfi'p'*)] 


dPn ' m ' ( cos  9) 

*  pnm^cos  0)  - ^0 -  m  COS  COS  m*^ 


dPnm ( cos  0) 


Pn ' m  * ( cos  6)®'  sin  mp  sin  m'$ 


We  thus  have 


r 

m  J  d<p  cos  m4  cos  m'$  «  m 5.,,  1 
0 


m'  /  d<f>  sin  sin  m'$  »  rnir6mm < 
0 


rrr  .  ♦(1.0)  +(2,e)  nm6nun'  b  (  (1)  .  a 

///  dT  ^pm  *  ^'p'm'  =  /a  dr  |r^np(3Ynp  r)  ^7 
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Hence  we  have 


///  a*  ^io1 

and  for  the  same  reasons  we  obtain  the 

///  ax  • 


•  -  o 


last  of  the  six  relations  we  need: 


+(2,o)  n 
En'p'm'  =  0  * 
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C  TEM 

COMMON  R.TEXT 

INTEGER  TEXT(2) ,T(2) ,T1 (2) 

DIMENSION  X(200) ,Y(200) ,RTAN(200)  ,X1 1 (10) ,X1F(10)  ,XCROSS(10) 
DIMENSION  I FLAG (15) 

DATA  IMARK/ '  *'/ 

10  WRITE  (1,20) 

20  FORMAT  (///'ENTER  TYPE  OF  RUN'/ 

- 'TM1 /TM2/TM3/TM4/TM5/TE1 /TE2/TE3/TE4 ' ) 

WRITE  (1  ,30)  T( 1 ) ,T(2) 

30  FORMAT  (’PREVIOUS  RUN  WAS:  ' , 2A2 /// 

-'(OR  ENTER  "0"  TO  QUIT)') 

READ  (1 ,40)  (T1 (I) ,1=1 ,2) 

40  FORMAT  (2A2) 

IF  (TI(I).EQ.'Q  ')  GO  TO  1000 
IF  (T1 (1 ) .NE. '  ')  GO  TO  50 

T1 (1 )=T(1 ) 

T1 (2)=T(2) 

50  T(1 )=T1 (1) 

T(2)=T1 (2) 

TEXT ( 1 )=T1 (1 ) 

TEXT (2)=T1 (2) 

IF  ((T(1) .EQ. 'TM') .AND. (T(2) .EQ. '1  '))  ICALL=1 
IF  ( (T( 1 ) .EQ . ' TM ' ) .AND. (T(2) .EQ . '2  '))  ICALL=2 
IF  ( (T( 1 ) .EQ. 'TM') .AND. (T(2) .EO. '3  '))  ICALL-3 
IF  ((T(1) .EQ.'TM') .AND.(T(2) .EQ. '4  '))  ICALL=4 
IF  ( (T( 1 ) .EQ . ' TM ' ) .AND. (T(2) ,EQ . '5  '))  ICALL=5 
IF  ((T(1).EQ.'TE’).AND.(T(2).EQ.'1  '))  ICALL=6 

IF  ((T(I).EQ.'TE') .AND. (T(2) .EQ. '2  '))  ICALL=7 
IF  ( (T( 1 ) .EQ . ' TE ' ) .AND. (T(2) .EQ . ' 3  '))  ICALL=8 
IF  ( (T( 1 ) .EQ . ' TE ' ) .AND. (T(2 ) \EQ . '4  '))  ICALL=9 
60  WRITE  (1 , 70)  R 

70  FORMAT  ('ENTER  R  (PREVIOUS  VALUE  WAS  \F7.5,'  )') 

READ  (1,80)  R1 
80  FORMAT  (F10.8) 

IF  ( R1  .GE.1  .)  GO  TO  60 
IF  ( R 1  . EQ . 0  . )  R1  =  R 

R-R1 

NPNTS=200 
DELX= . 1 
K=0 
L=0 
M=0 

KEND=0 

DO  90  1  =  1  ,15 
I FLAG ( I  )  =0 
90  CONTINUE 

100  IF  (L.NE.1)  GO  TO  110 
DO  220  M=1 , 10 

IF  (M.GE.10)  GO  TO  430 
110  IF  (L.NE.O)  DELX=(XF -XI ) /NPNTS 
DO  210  1=0  .NPNTS 
X ( I ) -FLOAT ( I ) / 10 . 

IF  (L.NE.O)  X(I)=XI+FLOAT(I)*DELX 
XT=X ( I ) 

CALL  FNCTN  ( ICALL, XT , YT) 

Y(I)=YT 
A=(1 . -R) *X ( I ) 

S=SIN (A) 

C»COS(A) 
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RTAN(I)«S/C 
I CROSS -0 
DIFF-YT-RTAN(I) 

IF  (DIFF.LE.O)  ND-0 
IF  (DIFF.GT.O)  ND-1 
IF  (I.EQ.O)  ND1-ND 
IF  (ND.NE.ND1 )  ICROSS-1 

IF ( ( RTAN (I) .LT.O.) .AND. (RTAN(I-1 ) .GT.O.))  ICROSS-O 
I F ( ( RTAN ( I ) .LT.O.) .AND. (RTAN(I-1 ) .GT.O.) .AND. (RTAN(I-1 ) .LT. 
-Y( I - 1 ) ) )  ICROSS-1 

IF  (I.EQ.O)  ICROSS-O 
ND1-ND 

IF  (L.NE.O)  GO  TO  120 

IF  (ICR0SS.EQ.1)  K-K+1 

IF ( (K .EQ . 1 ) .AND. (ABS(DIFF) .LT .0 .0000 1 ) )  K=0 
IF  (K.EQ.O)  ICROSS-O 
IF  (ICR0SS.EQ.1)  X1I(K)=X(I ) -DELX 
120  IF  (L.NE.1)  GO  TO  130 

IF  ( ( I CROSS .EQ . 1 ) .AND . (ABS (DIFF) .GT . 10 . ) )  IFLAG(K)=1 
130  IF  ((I.EQ.O) .AND. ( ICROSS .EQ . 1 ) )  IFLAG(K)=1 
IF  ((L.EQ.1) .AND. (IFLAG(K) .EQ.1))  GO  TO  150 
IF  (ICROSS. EQ.1)  WRITE  (1,140) 

140  FORMAT  ('*') 

IF  (Y(I)*RTAN(I) .LT.O.)  GO  TO  150 

IF  (L.EQ.2)  GO  TO  170 

IF  (ABS(DIFF) .LE. 0.00001)  X1F(K)=X(I) 

IF  (L .EQ .0)  GO  TO  150 
IF  (K.EQ.O)  GO  TO  150 

IF  (ABS(Y(I)-RTAN(I)) .LE.O. 00001)  GO  TO  320 
IF  (ABS(X(I)-X(I-I)) .LE.O. 000001)  GO  TO  320 
150  IF  (ICROSS. EQ.1)  XI I(K)-X(1)-DELX 

160  IF  ((L.EQ.1) .AND. (ICROSS. EQ. 1))  X1F(K)«X(I) 

IF  ((L.EQ.O) .AND. (ICROSS. EQ.1))  X1F(K)-X(I) 

IF  ((L.EQ.1) .AND. (ICROSS. EQ.1))GO  TO  170 
IF  (L.NE.O)  X1F(K)«X(I) 

170  IF  (L.EQ.1)  GO  TO  200 

IF  (ICROSS. NE.1)  WRITE  (1,180)  X ( I ) , Y( I) , RTAN ( I ) 

IF  (ICROSS. EQ.1)  WRITE  (1,190)  X(I) ,Y(I) , RTAN ( I ) 

180  FORMAT  (FI 2 .5 , 5X . 2F1 2 . 5) 

190  FORMAT  (FI 2 .5 , 5X , 2F1 2 . 5 , 5X , ' * ' ) 

200  IF  ((L.EQ.1) .AND. (ICROSS. EQ.1 ))  XI=X1I(K) 

IF  ((L.EQ.1) .AND. (ICROSS. EQ. 1))  XF-XIF(K) 

IF  ((L.EQ.1) .AND. (ICROSS. EQ.1))  GO  TO  100 
IF  (L.EQ.O)  KEND-K 
210  CONTINUE 

IF  (L.NE.1)  GO  TO  230 
220  CONTINUE 

230  IF  (L.EQ.1)  GO  TO  310 
IF  (L.EQ.2)  GO  TO  350 
240  NFLAG-0 

DO  250  K-1.KEND 

NFLAG-NFLAG+IFLAG(K) 

250  CONTINUE 

K1-KEND-NFLAG 
WRITE  (1 ,260)  T( 1 ) ,T(2) ,R 
260  FORMAT  (// 1 3X , 2A2 ,’ ROOTS  FOR  R-’,F7.5// 

-10X,  ' K'  ,  1  3X ,  ' XI '  ,  10X,  ’XF7) 

DO  290  K-1 , KEND 

IF  (IFLAG(K) .EQ.O)  WRITE  (1,270)  K,X1 I (K) ,X1F(K) 

IF  (IFLAG(K) .EQ.1)  WRITE  (1,280)  K,X1 I (K) ,X1F(K) 
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270  FORMAT  ( 10X , 11 , 5X , 2F1 2 . 5) 

280  FORMAT  (10X.I1  ,  5X  ,  2F12  .5 , 5X, ' # ' ) 

290  CONTINUE 

WRITE  (1 ,300)  K1 

300  F0RMAT(//’ THERE  ARE  ',12,  '  CROSSOVERS’/ 

-'ENTER  #  OF  CROSSOVER  TO  EXPAND'/ 

-'OR  RETURN  TO  CONTINUE') 

READ  (1,370)  K 
IF  (K.EQ.O)  GO  TO  350 
310  XI-X1 I (K) 

XF-X1 F (K) 

L-1 

M-0 

IFLAG(K)=0 
GO  TO  100 

320  WRITE  (1 ,330)  K 

330  FORMAT  (///II,'  CROSSOVER  OCCURS  AT') 

WRITE  (1,190)  X ( I ) ,Y(I) ,RTAN(I) 

WRITE  (1,340) 

340  FORMAT  (///) 

GO  TO  240 
350  WRITE  (1,360) 

360  FORMAT  (//'ENTER  1  TO  START  OVER'/ 

2  FOR  EXPANDED  SCALE’/ 

-'OTHERWISE  RETURN  FOR  GRAF') 

READ  (1,370)  IGRAF 
370  FORMAT  (13) 

L-1 

M-0 

IF  (IGRAF.EQ.1)  GO  TO  10 
IF  ( IGRAF. EQ. 2)  GO  TO  380 
CALL  GRAF  (NPNTS , X , Y , RTAN) 

GO  TO  240 
380  WRITE  (1 ,390) 

390  FORMAT  ('ENTER  K’) 

READ  (1,370)  K 
WRITE  (1,400) 

400  FORMAT  ('ENTER  INITIAL  AND  FINAL  VALUES  OF  X’) 
READ  (1,410)  XI, XF 
410  FORMAT  (2F9.5) 

IF  ((XI .NE.O.) .AND. (XF.NE.O.))  GO  TO  420 
XI-XII(K) 

XF-XIF(K) 

420  L-2 
M-0 

GO  TO  100 
430  WRITE  (1 ,440) 

440  FORMAT  ('NOT  A  REAL  CROSSOVER') 

M-0 

I F  LAG  (  K  )  - 1 
GO  TO  240 
1000  STOP 
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RSQ«R*R 

IF  (ICALL.EQ . 1 )  GO  TO  10 
IF  (ICALL.EQ. 2)  GO  TO  20 
IF  (ICALL.EQ. 3)  GO  TO  30 
IF  (ICALL.EQ. 4)  GO  TO  40 
IF  (ICALL.EQ. 5)  GO  TO  50 
IF  (ICALL.EQ. 6)  GO  TO  60 
IF  (ICALL.EQ. 7)  GO  TO  70 
IF  (ICALL.EQ. 8)  GO  TO  80 
IF  (ICALL.EQ. 9)  GO  TO  90 

10  TM1N-(1 ,-R)*XT*(1 .+R*XSQ) 

TM1D«RSQ*XT**4-(1 . -R+RSQ)*XSQ+1 . 

YT“TM1N/TM1 D 
GO  TO  100 

20  TM2N=3*(1 . -R)*XT*(RSQ*XSQ**2-2*( 1 .-R)**2*XSQ+1 2 .) 

TM2D«(R*XSQ)**3-3*R*(2*RSQ-3*R+2.)*XSQ**2-18*(RSQ-2*R+1 .)*XSQ 
-+36. 

YT-TM2N/TM2D 
GO  TO  100 

30  A»15.-2*XSQ 

B=1 5 . -2*XSQ*RSQ 

i  *y<;o+4S 

D=(XSQ*RSQ)**2-21 *XSQ*RSQ+45 . 

TM3N-3*XT*(A/C-R*B/D) 

TM3D-1 .+9*R*XSQ*A*B/(C*D) 

YT=TM3N/TM3D 
GO  TO  100 

40  A«420.-55*XSQ+XSQ**2 

B-420 . -55*XSQ*RSQ+(RSQ*XSQ) **2 

C*420 .-1 95*XSQ+1 0*XSQ**2 

D-420 . - 1 95*XSQ*RSQ+1 0* (XSQ*RSQ) **2 

TM4N-XT*(A/C-R*B/D) 

TM4D-1 ,+R*XSQ*(A*B) / (C*D) 

YT-TM4N/TM4D 
GO  TO  100 

50  A=*4725  .-630*XSQ+1  5*XSQ**2 

B»4725.-630*RSQ*XSQ+15*(RSQ*XSQ)**2 

C-4725  .-2205*XSQ+1 20*XSQ**2-XSQ**3 

D-4725 . -2205*RSQ*XSQ+1 20*(RSQ*XSQ)**2-( RSQ*XSQ)**3 

TM5N-XT*(A/C-R*B/D) 

TM5D-1 .+R*XSQ*A*B/(C*D) 

YT-TM5N/TM5D 
GO  TO  100 

60  TE1 N“( 1 . -R) *XT 

TE1D-1 .+R*XSQ 
YT-TE1N/TE1D 
GO  TO  100 

70  TE2N-(1 . -R)* ( 1 .+(1 ./3.)*R*XSQ)*XT 

TE2D-1 .-(1 .-3*R+RSQ)*XSQ/3+RSQ*XSQ*XSQ/9 

YT-TE2N/TE2D 

GO  TO  100 

80  A-15.-XSQ 

B-1 5 . -RSQ*XSQ 
C-5.-2*XSQ 
D-5.-2*RSQ*XSQ 
TE3N “XT* (A/C-R*B/D) 

TE3D-3 ,+R*XSQ*A*B/ (3*C*D) 

YT-TE3N/TE3D 
GO  TO  100 

90  Y«XT*R 
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Y2-Y*Y 
Y3-Y2*Y 
Y4-Y3*Y 
Y5-Y4*Y 
X2-XSQ 
X3-X2*XT 
•X4=X3*XT 
X5-X4*XT 
A=105*XT-10*X3 
B=105-45*Y2+Y4 
C=105-45*X2+X4 
D-105*Y-10*Y3 
TE4N=A*B-C*D 
TE4D=C*B+A*D 
YT-TE4N/TE4D 
100  RETURN 
END 


SUBROUTINE  GRAF  (NPNTS , X , Y1 , Y2) 
COMMON  R.TEXT 
INTEGER  TEXT (2) 

INTEGER  RTEXT(3) 

DIMENSION  X(200) ,Y1 (200) ,Y2(200) 
ENCODE(6 , 5 , RTEXT)R 
5  FORMAT  ('R-' ,F3.2) 

CALL  MINMAX  (Y1 , NPNTS, 1 , Y1MIN , Y1 MAX) 
IF  (Y1MAX.GT.150.)  Y1MAX-150. 

IF  (Y1MIN.LT. -99.)  Y1MIN--99. 

DO  10  1-1 .NPNTS 
X(I)-X(I) 

IF  (Y1 (I) .GT. 150  .)  Y1 ( I)-1 50  . 

IF  (Y1  (I)  .LT.-99  .)  Y1(I)  — 99. 

10  CONTINUE 
CALL  SETPDQ 
CALL  SCREEN 
CALL  IN IT  (1.5,1  .5) 

CALL  SCALE  (X, NPNTS ,8 . , 1 ,0) 

CALL  SCALE  (Y1 .NPNTS , 5 . , 1 , 1 ) 

CALL  ENTGRA 

CALL  XAXIS  CX’  ,1,8.) 

CALL  YAXIS  ('FUNCTIONS  FOR  R',15,5.) 
CALL  DATAQ  (X ,Y1 .NPNTS , 1 , 1 ) 

L-3 

DO  20  1-1 .NPNTS 

CALL  DRAWC  (X( I ) , Y2 ( I) , L, 4) 

L-4 

20  CONTINUE 

CALL  GRID  (8  . , 5 . , 1 . , 1 . ,4 ,0) 

CALL  ABSVEC  (1  .7,6.  ,0) 

CALL  SYMBOQ  (6 .RTEXT , 7) 

CALL  ABSVEC  (8. 6, 6. 7,0) 

CALL  SYMBOQ  (4, TEXT, 7) 

CALL  EXITGR 
CALL  WAIT 
CALL  SCREEN 
RETURN 
END 
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